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Summary 


A discrete-vortex method has been developed to model the vortex flow over a delta wing and to 
investigate the wing rock problem associated with highly swept wings. The method uses two logarithmic 
vortices placed above the wing to represent the vortex flow field, but instead of the boundary conditions 
used in previous discrete-vortex models, it uses conditions based on conical flow, vortex rate of change of 
momentum, and other considerations to position the vortices and determine their strengths. For symmet- 
ric zero roll, the assumptions are that the velocity at the center of each vortex is directed radially outward 
from the wing axis (a conical flow requirement), that the force on the vortex system is equal to the rate of 
change of momentum of the vortex system, and that the maximum force on the vortex is equal to that 
possible with attached flow. For asymmetric nonzero roll, the strengths and positions of the vortices were 
determined by using additional assumptions about positions of the vortices relative to streamlines and 
lines of constant velocity potential. A strictly kinematic relationship based on the time analogy and 
conical-flow assumptions is derived in an appendix to this report and was used to determine the hysteretic 
positions of the vortices during rolling oscillations. The rolling moments calculated for the wing and 
vortex system were used in the single-degree-of-freedom equation of motion for the wing about its roll 
axis to generate wing rock time histories. The calculated results in this report are compared with avail- 
able experimental data and analyzed to show the effects of angle of attack, wing geometry, and wing 
inertia on wing rock characteristics. 

The results show that the method was able to model the basic features of wing rock once the static 
rolling moment characteristics were adjusted to agree with the experimental data. Static and dynamic 
vortex positions and wing rock amplitudes and frequencies were generally in good agreement with the 
available experimental data. The results verify that wing rock is caused by hysteretic deflections of the 
vortices and indicate that the stabilizing rolling moments that limit the amplitude of wing rock oscillations 
are essentially the result of one primary vortex moving outboard of the wing to where it has less influence 
on the wing. 

Introduction 

One of the more serious problems affecting fighter aircraft at high angles of attack is the condition 
known as wing rock — a high amplitude, high frequency rolling oscillation associated with the leading- 
edge vortex flow on highly swept wings. Wing rock can cause serious maneuvering and tracking prob- 
lems and, because it involves a substantial loss of lift, also can present a potential safety hazard during 
takeoff and landing for both fighter and highly swept transport designs. 

Although wing rock has been studied in a number of investigations (see, for example, refs. 1-16), it is 
still not completely understood. The experimental investigations (refs. 1-9) have indicated that wing rock 
is probably caused by the time lags in vortex position that occur when a wing with leading-edge vortex 
flow is given a roll velocity. As the wing rolls, the primary vortex on one side of the wing will move 
closer to the wing upper surface, while that on the other side moves away from the wing upper surface, 
producing an unstable rolling moment that drives the motion to higher amplitudes. The experiments, 
however, have not been able to show what causes the restoring moments that come into play at the higher 
roll angles and limit the amplitudes of wing rock oscillations. Also not known is how parameters such as 
wing size, sweep, and moment of inertia affect wing rock amplitudes and frequencies. 

Computational studies of wing rock have not added much information beyond what has been learned 
experimentally. Navier-Stokes calculations for the wing rock problem are extremely difficult since they 
involve unsteady motions of both the wing and flow field. Although Navier-Stokes calculations of wing 



rock were made in reference 1 0, they were for the forced oscillation of a delta wing at an angle of attack 
of 20°, an angle for which wing rock would not be expected to occur (see, for example, refs. 4, 7, and 9). 
Two-dimensional Euler calculations for a delta wing undergoing wing rock at supersonic speeds were 
made in reference 11, and fully three-dimensional Euler calculations were made in reference 12. These 
Euler calculations were able to simulate limit-cycle wing rock oscillations and show that they were the 
result of moment variations that were unstable at the lower wing roll angles and stable at the higher roll 
angles; however, neither was able to identify the underlying aerodynamic mechanisms that cause this type 
of rolling moment variation. Moreover, neither the Euler nor the Navier-Stokes calculations made to date 
have produced results that can be directly compared with available experimental data. 

The discrete-vortex models developed in references 1 3 and 1 4 represent a much simpler method for 
studying wing rock. These models assume that the vortex flow field is dominated by two primary vor- 
tices that form above the wing and which can be represented mathematically by simple logarithmic func- 
tions. They neglect the secondary and tertiary vortices that fonn near the wing leading edges in a real 
flow and assume that vorticity is fed into the primary vortices along straight-line feeding sheets connect- 
ing the wing leading edges to the vortex centers, rather than the curved feeding sheets assumed in the 
more complex discrete-vortex models. Besides greatly simplifying the problem, this type of model pro- 
vides better insight into the aerodynamic mechanisms involved in wing rock than the more complex 
methods. These models reduce the problem to an essentially three-body vibrations problem: once the 
positions and strengths of the vortices have been determined, the moments they produce on the wing can 
be determined and easily analyzed to provide a clearer picture of the aerodynamic mechanisms involved. 

The problem with discrete-vortex methods, at least in the simplified fonn in which they are developed 
in references 13 and 14, is that they are inherently inaccurate. They are based on the Brown and Michael 
model described in reference 16, which does not give accurate vortex positions when applied to two- 
dimensional cylinders and cones (ref. 17) or delta wings (ref. 16), even for the basic case in which the 
vortices are symmetrically aligned. The models developed for wing rock, therefore, can at best be 
expected to provide qualitative but not quantitative results. The inaccuracies in these discrete-vortex 
methods are usually thought to result from the simplifying assumptions made in their development; 
however, the problem is not that these models oversimplified the flow field, rather the problem is that the 
proper boundary conditions were not applied in obtaining a solution. In reference 1 8, it was shown that a 
more accurate discrete-vortex method could be developed for cones by using boundary conditions based 
on conical vortex flow and the momentum of the vortex system, rather than the conditions used in the 
basic Brown and Michael model. 

The purpose of the present paper is to show that by using a different set of boundary conditions than 
those used in previous discrete-vortex models, along with other constraints on vortex motion, it is possi- 
ble to develop a more accurate method for studying wing rock. This new method not only can provide a 
more thorough understanding of the wing rock problem, but also can give reasonably accurate predictions 
of wing rock amplitudes and frequencies. 

The method developed in this paper discards the two basic boundary conditions used in the Brown and 
Michael model (ref. 16) and the wing rock models of references 13 and 14. One is a Kutta condition, 
which requires the flow to leave smoothly from the wing leading edges; this condition is considered inap- 
plicable because the present method neglects secondary and tertiary vortex flow. The second is a zero- 
force condition in which the force on an assumed feeding sheet is set equal to the force caused by induced 
effects at the vortex center. It replaces these two conditions with a series of others designed to determine 
the vortex positions and strengths during different stages of the wing rock process. 
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In the present method, the problem is solved first for symmetric zero roll; this is done by applying a 
conical flow condition, which sets the strengths of the vortices, and then applying a momentum condition, 
which equates the force on the vortex system to its rate of change of momentum. A maximum force con- 
dition that equates the force on the vortex system to the maximum force possible with attached flow is 
then applied as the third condition needed to solve the zero roll case. 

Information about the streamlines on which the primary vortices lie in the symmetric case is then used 
to determine vortex positions and strengths with the wing at a nonzero roll angle. The vortices are 
assumed to move with these streamlines when the wing is rolled and also to retain their positions relative 
to certain velocity-potential lines that remain fixed with respect to the wing. These conditions, along with 
the conical flow and momentum conditions, provide the six conditions needed to solve the problem for 
asymmetric nonzero roll. 

The final positions that the vortices assume when the wing undergoes a rolling velocity are determined 
by adding incremental hysteretic displacements to the calculated static positions. An equation that is 
based on the conical flow and time analogy assumptions and gives the hysteretic displacements in terms 
of static vortex position, wing roll rate, semiapex angle, and angle of attack is derived in this report. 

Once the positions and strengths of the vortices have been determined in the present method, then 
methods based on Blasius integrations (refs. 19 and 20) are used to determine the static and dynamic 
rolling moments acting on the wing. These are then used in the single-degree-of freedom equation of 
motion for the wing about its roll axis to generate wing rock oscillations. To assess the accuracy of the 
method, the computed results are compared with available experimental data for ranges of angle of attack, 
free-stream velocity, wing geometry, and moments of inertia. 


Symbols 


a 

b 

C F, 

C k 

Ci 

Cl 

c, 


1 attached 


b~ 

wing aspect ratio, — 

s 

local semispan of wing, ft 
wing span, ft 

. , ~ ,. ..... . Attached-flow force 

attached-tlow force coefficient, , lb 


qS 


nondimensional vortex strength, 
rolling moment coefficient. 


InaUco sin a 
Rolling moment 


attached-flow lift coefficient. 


qSb 

Attached-flow lift 

qS 


rolling moment coefficient due to friction in roll-oscillation apparatus used in wind tunnel 
test 
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Cl 


h 


rolling moment coefficient due to hysteretic effects, 


qSb 


Ci 


'j 

attached-flow roll damping coefficient. — — , 

( A*U \ 

a 


\2C/oo / 


C; static rolling moment coefficient, — — 

5 qSb 

Ci I static rolling moment coefficient calculated with vortex at static position Xk,Vk at given 

' x k*yk 

wing roll angle 


Ci 


\x k +Ax k ,y k 


static rolling moment coefficient calculated with vortex at hysteretic position x\ + A x^yk at 
given wing roll angle 


static spring rolling moment derivative, 


dQ 

rfcj) 


C N 

Co 

Cx 

Cy 


. r .... . Normal force 

wing nonnal force coefficient, 

qS 

path of integration around closed curve in /.-plane; positive direction is that for which an 
observer, traveling along path, would keep enclosed area to left 

X-axis force coefficient, 

qS 


7-axis force coefficient, 


qS 


c wing root chord, ft 

calc. calculated 

cl l, (h centroids for vortex pairs 

exp. experiment 

F steady pressure force on wing, lb 

Fps force on feeding sheet assumed in Brown and Michael model of reference 16, lb 

F max maximum force in y-direction on either wing or primary vortices (see fig. 6(d)), lb 

Fy force due to nonzero velocity at primary vortex center, lb 

Fx steady pressure force on wing in x-direction, lb 
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1 unsteady 


Fy 


1 unsteady 


/ 

g 


h 


h 

i 

k 


L 

l 

l h 

Ip 

Is 

MV 

P 


<7 


r C k 


r a 

S 

t 

U 


unsteady force on wing cross section in ^-direction, lb 

steady pressure force on wing in y-direction, lb 
unsteady force on wing cross section in y-direction, lb 

oscillation frequency, cycles/sec 

influence coefficient, which, when multiplied by vortex circulation for symmetric vortex 
alignment, gives X-axis velocity induced at given primary vortex center by other vortices 
in system 

influence coefficient, which, when multiplied by vortex circulation for symmetric vortex 
alignment, gives 7-axis velocity induced at given primary vortex center by other vortices in 
system 

wing moment of inertia about longitudinal (roll) axis, slug-ft- 
unit vector along imaginary axis 

constant used to account for difference between experimental rolling moments and calcu- 
lated, attached-flow rolling moments 

length of wing root chord, ft 

distance along wing root chord, positive from nose rearward, ft 

rolling moment caused by hysteretic effects, ft-lb 

rolling moment caused by wing rolling velocity, ft-lb 

static wing rolling moment, ft-lb 

linear momentum, lb-sec 

static pressure, lb/ft 2 

1 2 

free-stream dynamic pressure, — pU x , lb/ft 2 

I 9 1 

scalar distance to primary vortex center in a-plane, x +y , ft 

V c k c k 

radius of cylinder in a-plane, ft 
wing area, ft 2 
time, sec 

X-axis component of velocity that results from free-stream flow around wing at given cross 
section, nondimensionalized with respect to U, rj sin a 

free-stream velocity, ft/sec 
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u X-axis component of velocity at center of primary vortex, nondimensionalized with respect 

to Uao sin a 

V 7-axis component of velocity that results from free-stream flow around wing at given cross 
section, nondimensionalized with respect to U, 0 sin a 

V max maximum velocity along either cylinder or M* = 0 streamline in figure 6, ft/sec 

v 7-axis component of velocity at center of primary vortex, nondimensionalized with 

respect to U !tJ sin a 

W complex potential of flow field, <f» + /'l', ft2/sec 

W complex conjugate of W 

X coordinate axis normal to wing cross section (see fig. 1) 

x X-axis coordinate (see fig. 1 ) 

x c coordinate along real axis in circle plane (see fig. 2(a)) 

x Cf x c distance to centroid of vortex pair consisting of vortex k and its image (see fig. 9(b)) 

xj state variable, c|> 

x o state variable, — 

dr 

7 coordinate axis that lies in and is parallel to plane of wing cross section (fig. 1) 

y 7-axis coordinate (fig. 1) 

y c coordinate along imaginary axis in circle plane, fri/sec (see fig. 2(a)) 

a wing angle of attack, deg 

wing sideslip angle, deg 

I* circulation of vortex k, positive clockwise in first ( +x, +y) quadrant ofX-7 plane, ft 2/sec 

Y rolling moment reduction factor, defined by equation (42) 

Ax hysteretic deflection of vortex in x-direction, ft 

Ax/, hysteretic deflection at wing trailing edge of vortex in x-direction, ft 

8 wing semiapex angle, deg 

t, complex coordinate of primary vortex center in ^-plane, + ir\, ft 

r) coordinate along imaginary axis in 'C-plane, ft 
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0 


angle between real axis and line extending from axis origin through primary vortex center 
in a-plane, deg 


A wing sweep, deg 

X complex coordinate of primary vortex center in physical plane, x + iy, ft 

K s complex coordinate of separation point at wing leading edge in Brown and Michael model 

of reference 1 6, ft 

A complex conjugate o f a, ft 

5 coordinate along real axis in C-plane, ft 

p fluid density, lb-sec 2 /ft 4 

a complex coordinate of primary vortex center in circle plane, x c + iy c , ft 

a complex conjugate of a 

<t> velocity potential, equal to real part of complex potential, ft 2 /sec 

<t> ( j ) = 0 o velocity potential with wing at zero roll angle, IWsec 

<h^ =0 , velocity potential for vortex k with wing at zero roll angle, ft 2 /sec 

4) wing roll angle, deg 

4>max maximum amplitude of wing rock oscillation, deg 

'P stream function, equal to imaginary part of complex potential, ft 2 /sec 

W ^ =Q0 stream function defining streamline on which vortex k lies with wing at zero roll angle, 

ft 2 /sec 

Subscripts: 

k primary vortex index: 1 when referring to conditions at center of vortex 1; 2 when referring 

to conditions at center of vortex 2 (fig. 2) 

te wing trailing edge 

Mathematical Development 

Basic Assumptions 

The present method makes four basic assumptions about the flow field. First, the flow is assumed 
incompressible and inviscid. Second, all of the vorticity generated by flow separation at the leading edge 
of a delta wing at an angle of attack is assumed to be concentrated in two primary vortices that form 
above the wing and for which the viscous core region is small enough to be neglected. Third, the flow is 
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assumed to be conical. This third assumption allows the problem to be treated as a two-dimensional 
problem in the cross-flow plane of the wing (the axis system for the wing and cross-flow plane is shown 
in fig. 1 ) and plays an important role in establishing some of the boundary conditions needed in the solu- 
tion. The fourth assumption is that the two-dimensional time analogy applies. According to this analogy, 
changes that take place with distance, in going from cross section to cross section along the axis of a 
three-dimensional delta wing, are equivalent to those that would take place with time for a plate whose 
span is increasing with time in a two-dimensional flow field. (For an example of the use of the time anal- 
ogy, see ref. 21.) This analogy allows time in the unsteady two-dimensional flow to be related to distance 
in the three-dimensional flow as 


t/ooCosa 

In addition to these basic assumptions, other assumptions have been made throughout this report in 
connection with the development of boundary conditions and the determination of wing moments. These 
assumptions are indicated and discussed as they arise in the model development that follows. 


x 



Figure 1. Wing axis system. 


Basic Flow Model 

General Description 

Under the assumptions of an incompressible inviscid fluid, the governing equation for the flow is 
Laplace’s equation. The problem is to find a solution in the form of a complex potential that satisfies 
Laplace’s equation, subject to the appropriate boundary conditions. In the present method, the flow in the 
crossplane is first modeled in a complex cr-plane (fig. 2(a)), where the complex potential can be con- 
structed from known elementary potential functions and a tangent flow boundary condition can be applied 
to the body. The flow in the cr-plane is then transformed into flow about the wing in the physical /.-plane 
(fig. 2(b)) where the boundary conditions on the vortices are applied. In the cr-plane, the flow model con- 
sists of two logarithmic vortices, of strengths L| and Y 2 , placed behind a two-dimensional circular cylin- 
der (the cross-flow velocity in fig. 2 is from left to right). Image vortices, with strengths equal and oppo- 
site to those of the primary vortices, are placed inside the cylinder to satisfy the tangent-flow boundary 
condition at the surface. 
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(a) a-plane. 


ty 



direction 

(b) k -plane. 


Figure 2. Vortex system in a- and k-planes. 


Complex Potential for Wing-Vortex System 

For the most general case in which the vortices are asymmetrically aligned and the wing is rolled 
through an angle <jr, the complex potential for the flow about the wing-vortex system in the a-plane is 

w(o,<k) = <s>+m> 


= alien sina 


+ i C] In 


— (cos <() + / sin<()) - 


a 


(cos 4> + / sin cp) 


17 

\ 


( 

2 \ 


/ 

\ 


t 

2 Ml 

a Oj 


O 

_fo_ 

+ i C2 In 

o__o 2 


Q 

1 

'h" 


r o J 





\ r o 

r o y 



a ^JJ 


( 2 ) 


The first term in the first set of brackets on the right-hand side of equation (2) is the complex potential 
for a free-stream flow inclined at angle $ with respect to the real Tf c -axis. The second term in these 
brackets is the complex potential for a source-sink doublet that is also aligned at angle (j) with the real 
axis. The remaining tenns on the right-hand side represent the contributions of the primary vortices and 
their images. The complex potential given by equation (2) has branch points at the centers of each pri- 
mary vortex and its image and is, therefore, multivalued. It can be made single-valued — allowing a 
unique solution to be obtained for a given set of boundary conditions — by drawing branch cuts between 
the primary vortices and their images. 
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Coordinate and Velocity Transformations 


The flow in the physical /.-plane is for a flat plate, which at zero roll is aligned perpendicularly to the 
free stream. The equation for transforming coordinates from the a- to the /.-plane is 


X _ l(a _ 1 ' 
a 2 \r 0 a/r ()/ 


( 3 ) 


Transformations for going back and forth between any of the three complex planes used in the present 
analysis are given in appendix A. 

Flow velocities in the physical plane are obtained by taking the derivative of the complex potential 
with respect to X (the coordinate in the physical plane). This derivative can be written as 


U-iV- 


1 


dW 


1 


dW da 


t/ooSina dX U^sinci do dX 


(cosA + isind)') 

V 7 cr 


(cos 4> + / sin 4>) 


+ i C] 


O Or 


a r‘ 


\ r o at ) 


+ i C't 


a a 2 


a r: 


\ r o °2 ) 


1 + 


X 

a 


X z 


+ 1 


( 4 ) 


Note that when equation (4) is used to obtain the velocity at a given vortex center, the velocities 
induced by that vortex on itself are neglected. For example, when the velocity at a = o\ (the center of 
vortex 1) is to be determined, the tenn 


1 


a a j 

r, K, 


( 5 ) 


' o ' o 

which otherwise would be infinite, is set equal to zero. 


Solution Procedure 


The equation for the complex potential (eq. (2)) contains six unknowns: x\, y\, x 2 , >'2 (the coordinates 
of the centers of the two primary vortices), C\, and C 2 (the nondimensional strengths of the vortices). 
These positions and strengths can be determined by applying the appropriate boundary conditions; once 
they have been determined, the rolling moments caused by the vortices can be computed and used in 
the equation of motion for the wing-vortex system to generate wing rock oscillations. In the present 
method, the static vortex positions and strengths are determined first for symmetric zero roll by using 
motion and force restraints on the vortices, along with the equation relating the force on the vortex system 
to its rate of change of momentum. Information about the streamlines on which the vortices are located 
for the symmetric roll angle is then used to obtain vortex positions and strengths for the asymmetric 
nonzero roll angles. Flowcharts outlining the procedures for determining the vortex positions and 
strengths for zero and nonzero roll angles are presented in figures 3(a) and (b), respectively. These can be 
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used as the basis for generating computer codes for performing the calculations required in the present 
method. The boundary conditions needed to apply these procedures are stated and discussed in the fol- 
lowing sections. 

Determination of Static Vortex Positions and Strengths at Zero Roll Angle 

With the wing at zero roll angle, the vortices are symmetrically aligned about theX-axes, so that 


*2 = x l 


y 2 = -3T ■ 
C 2 =-Cj 


( 6 ) 


Therefore, only the coordinates and circulation of one of the primary vortices must be determined. In the 
discrete -vortex methods of references 13 and 14, these three unknowns were determined by imposing a 
zero-force condition (which provided two of the needed equations) and a Kutta condition (which provided 
the third). In the present investigation, these unknowns were determined by applying three symmetric 
boundary conditions (S-l, S-2, and S-3) involving vortex motions and forces. 


Symmetric Boundary Conditions for Wing at Zero Roll Angle 

S-l. The velocity at the primary vortex centers is in the radial direction. This condition allows non- 
dimensional vortex strength to be determined according to the equation: 

Cl .-C,.-M« (7) 

h \-g\(y i /■*]) 


S-2. The force on the vortex system in the x-direction (normal to the wing) is equal to the rate of 
change of momentum of the vortex system in the x-direction, a condition expressed by the equation: 




• 

' 1 dW ^ 

f /M 


Re 

i ( 

> 


\d - 

= 4jtCi 




v (Zoo sina dk j 

1 \a ) 



r Cl > r o j 


sin 0j 


tan 6 
tana 


( 8 ) 


where the integration in equation (8) is carried out around a path C 0 that includes only the two primary 
vortices. 


S-3. The maximum lateral force for the flow is equal to the maximum lateral force that would be 
exerted at the wingtip in attached flow, a condition expressed by the equation: 


i . 



( 1 dW ^ 

\ 2 Jk\ 

- Im 

c 

) 


\ d — 

2 


r 

^C/ooSina cfk j 

\o / 


( 9 ) 


The term on the left-hand side of equation (9) is the absolute value of the computed force coefficient; the 
path of integration is taken around only one of the primary vortices and not around any other singularity 
in the flow field. 

These three boundary conditions are applied according to the procedure outlined in figure 3(a), and are 
discussed in detail in the next sections. 
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(a) Zero roll angle. 

Figure 3. Flowcharts of procedure for determining static vortex positions and strengths. 
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Discussion of Symmetric Boundary Conditions 

Radial-Velocity Condition (S-l). The radial-velocity condition is one of the conditions imposed for 
cones in reference 1 8 and determines the vortex circulation for a given vortex location. It is based on the 
fact that the conical flow assumption for a delta wing requires the vortex centers to move radially away 
from the centerline of the three-dimensional wing as the vortices progress from one cross section of the 
wing to the next. Since the primary vortices are assumed to be free vortices, each moving in the direction 
of the local flow velocity at its center, this condition is satisfied if the local induced velocity at the center 
of each primary vortex is directed radially outward, along a line that extends from the center of the plate 
through the vortex center (see fig. 2(b)). 

To determine the vortex circulation using this condition, the x- and y-components of the velocity (see 
eq. (4)) at the center of one of the primary vortices are first written in the form of a free-stream compo- 
nent plus an influence coefficient times the nondimensional circulation of the vortices as follows: 


v l =V ] +h l C l 

(TO) 

u \ = U\ + g\C\ 

(11) 


Then if the velocity is radial, the slope v\/u\ of the velocity vector must equal the slope y\!x\ of the line 
running from the origin through the vortex center; that is: 


Collecting terms gives 


v’l V i +h ] C ] >’! 

«1 U\+g\C\ x \ 


h \~S\ 


Zi 
x \ ) 


Cl =- 


Vi-Ui 


Zi 
X 1 ) 


Solving for Ci gives equation (7). 


( 12 ) 


(13) 


Applying the radial-velocity condition will give the nondimensional vortex strength for any assumed 
value of .Yi and vi- This operation is perfonned in the top part (down through the element Ci) of the flow- 
chart in figure 3(a). Contour lines showing the positions that the vortices would take for specified values 
of nondimensional vortex strength are shown in figure 4. 


A-Axis Momentum Condition (S-2). Equation (8) is the equation of motion for the vortex system. It 
replaces the zero-force condition of the Brown and Michael (ref. 16) and Bryson (ref. 17) models in 
which it was assumed that the force on the straight-line feeding sheet connecting each primary vortex to a 
separation point on the body is balanced by the force caused by induced effects at the vortex center. In 
the Brown and Michael and related models, these feeding sheets serve as the branch cuts needed to make 
the complex potential for the vortex system single-valued, and, as such, are subjected to a force (see 
ref. 17). Across such a branch cut, the velocity potential <I> changes by the amount T, the vortex circula- 
tion, so that when T changes with time, the pressure difference across the sheet is 


, d<l> 
dp = P— = P 
dt 


dr k 

dt 


(14) 
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Figure 4. Lines of constant circulation in X-Y plane. 

When integrated over the length of the feeding sheet, this pressure difference results in a pressure force of 

F FS =- ip ^t(X k -X s ) (15) 

In the Brown and Michael model, this pressure force is balanced by a force resulting from a nonzero 
velocity at the vortex center, a force given by the Kutta-Joukowski theorem as 

^ = <pr*(v,-^j (16) 

Setting the sum of forces Fpg and Fy equal to zero gives the zero-force conditon used in the basic Brown 
and Michael model as 

-,pi&p t -x s ) +I pr i (v i -^)=o 07, 

Although the Brown and Michael zero-force condition has been used extensively in discrete-vortex 
models, it is not a valid condition. As indicated in appendix B, the force on the feeding sheet is an un- 
steady force caused by changes in vortex strength and, as such, should not be balanced by a steady force 
at the vortex center. 

Instead of using the Brown and Michael condition, the present method specifies a force balance by 
equating the X-axis force on the vortex system to the rate of change of momentum in the x-direction of the 
vortex system. (In the v-direction, the net force and the rate of change of momentum of the vortex system 
are both zero for the zero roll case; therefore a 7-axis momentum equation would not provide any useable 
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information.) In the x-direction, the force on the vortices obtained by Blasius integration (see eq. (Bl) of 
appendix B) can be expressed in nondimensional fonn as 


Re 



' 1 dW 

^I/ooSina cTk 



(18) 


where the integration is carried out around a closed path C a , taken in the counterclockwise direction, that 
encircles the two primary vortices but not any other singularities. The force on the wing, which can be 
obtained by integrating around a path that includes the wing but not the vortices, will be equal and 
opposite to this force on the vortices. 

The rate of change of momentum of the vortex system is determined in the a-plane (see fig. 2(a)), 
where the primary vortices appear as pure logarithmic vortices whose complex potential has the fonn 

W(a) = ^\n(o-o k ) (19) 

2ji 


For a pair of vortices with this mathematical fonn, the momentum is equal to pT times the distance 
between the two vortices. That is, in the a-plane, the momentum in the direction of the real axis for a 
vortex pair consisting of the upper primary vortex and its image is equal to 


/ 

MV = ipT* 

V 


.2 \ 


r k~- 


r k 


sin 0/. 


( 20 ) 


In the direction of the real axis, the total momentum for the two pairs forming the vortex system is twice 
the momentum calculated in equation (20). 


The rate of change of momentum of the vortex system can be written in terms of conditions at 
vortex 1, for which the circulation is positive, as 


d(MV ) 
dt 




2 \ 




n n, dT 1 

sin 0i + 2p — - 
dt 



r > ) 


sin 9] 


( 21 ) 


As discussed in appendix B, the unsteady force caused by changes in vortex circulation (the second 
term on the right-hand side of eq. (21)) is exactly equal to the rate of change of momentum caused by 
changes in vortex circulation, regardless of where the vortex center is located. No new information can 
be obtained, therefore, by equating force caused by changes in circulation to the rate of change of 
momentum caused by changes in circulation. In the present analysis, therefore, neither the unsteady force 
caused by changes in vortex circulation nor the change in momentum resulting from changes in circula- 
tion is considered when determining vortex position. 


Equation (21), without the tenn for rate of change of circulation, can also be written as 


d(MV) 

dt 


= 2p Tj 




sin 0 


dr Q da dl 
da dl dt 


( 22 ) 
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By making the substitutions r Q = u/2, da/dl = tan 8, and dl/dt = U x cos a, equation (22) can be written as 


d(MV) 

dt 


= P r i 



. 2 \ 


sin 0]tan6 U^cosa 


(23) 


In evaluating the derivative with respect to r 0 in equation (23), it should be noted that r\ can be written 
as r\ = r 0 (r\/r 0 ), and that the quantity (r\/r 0 ) is constant for each station of the delta wing in conical flow; 
hence, dr\!dr 0 = r\!r 0 . Thus, equation (23) can be written as 


d(MV) 

dt 



I S 

r \l r o, 


sin 0] tan 5 U x cos a 


In nondimensional fonn this becomes 


1 d(MV) 
p Ula dt 



1 N 

r \l r o, 


Cj sin 0] 


tan8 . ? 

sm" a 

tana 


(24) 


(25) 


Equating this nondimensional change in momentum to the force coefficient given by equation (18) yields 
the momentum equation for the symmetric vortex system, given by equation (8). 

For a given value of tan a/tan 6, equation (8) will be satisfied for certain values of x\ and y\, the vor- 
tex coordinates in the physical plane. As shown in the flowchart of figure 3(a), once initial values of x\ 
and _v'i have been specified, C\ can be determined from the radial-velocity boundary condition (S-l), and 
x\ can be iterated to get the final values of jq and y\ that will satisfy equation (8). Contours of constant 
tan a/tan 8, obtained by applying the X-axis momentum condition, are shown along with the constant Q- 
contours, obtained by applying the radial-velocity condition, in figure 5. 

Maximum Lateral Force Condition (S-3). Conditions S-l and S-2 would provide the information 
needed to determine vortex positions and strength — if the circulation produced by the flow separation at 
the wing leading edge were known. For cones in reference 18, this circulation was calculated as being 
proportional to the square of the edge velocity at the separation points on the cone, and although the loca- 
tions of the separation points on the cone surface could not be determined using the method of refer- 
ence 18, they could be approximated with available experimental data. The velocities at these separation 
points, as well as at all points on the cone surface, are finite, and the circulation being made available for 
vortex formation could be calculated. 


For delta wings, on the other hand, although the flow is known to separate at the sharp leading edges 
of the wing, the velocity at the separation points, at least in the present model, is infinite. Therefore, sepa- 
ration point velocities cannot be used to calculate the vortex circulation. Condition S-3 is a way of 
working around this problem. Instead of determining vortex positions and strengths by calculating the 
circulation made available for vortex formation, they are determined by putting a limit on the maximum 
force that can be produced on the wing by the flow field around it. This condition assumes this force 
cannot exceed the leading-edge suction force that the wing would develop with completely attached flow. 

Condition S-3 is strictly an assumption, presented without proof, but is analogous to the maximum- 
velocity condition imposed for a two-dimensional circular cylinder with trailing vortex flow in 
reference 22. Figure 6 illustrates this analogy by presenting attached- and separated-flow velocity 
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Figure 5. Lines of constant circulation and of constant tan a/tan 8. 


distributions for the circular cylinder of reference 22 (figs. 6(a) and (b)) with typical attached- and 
separated-flow force distributions computed for the delta wing cross section by Blasius integrations 
(figs. 6(c) and (d)). 

To obtain solutions in reference 22, the assumption is that when the flow behind the cylinder separates 
and forms into vortices, the maximum flow velocity possible with vortex flow should be the same as that 
for the cylinder with attached flow. As shown in figure 6(a), with attached flow, the flow velocity 
reaches a maximum value of 2 f/ x . at the points x = 0, y/r Q = ±1 on the surface of the cylinder; with sepa- 
rated flow (fig. 6(b)), it reaches a maximum at a point that lies above and slightly behind the vortex center 
on the 'F = 0 streamline that surrounds each vortex. In reference 22 the assumption is that this maximum 
velocity, regardless of where it occurs in the vortex flow field, must still be equal to 2U<n as shown in 
figure 6(b). For this velocity to be higher than 2 would mean that the vortices were adding energy to 
the outer flow instead of deriving their energy from it. 

Condition S-3 is a similar condition but places a limit on the forces acting at the delta wing cross sec- 
tion, rather than on the velocities. The maximum vortex force for a delta wing with attached flow is that 
given by the Polhamus leading-edge suction analogy (ref. 23). This force can be determined by Blasius 
integrations around half the wing cross section, and for a flat-plate cross section it has a value at each 
leading edge (in coefficient form) of ( nil) sin 2 a, as shown in figure 6(c). (These forces, according to the 
Polhamus leading-edge suction analogy, turn 90° to become the normal force on the wing.) This force 
arises from the fact that although there are infinite velocities and pressures at the leading edges of the 
wing with attached flow, they act over an essentially zero area to produce a finite lateral force. 

With separated flow (fig. 6(d)), there is still an infinite velocity but a finite force at the leading edge 
because a Kutta condition is not being satisfied in the present model; however, these calculated forces are 
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(c) Delta wing cross section with attached (d) Delta wing cross section with separated vortex 

flow; a = 30°. flow; a = 30°. 

Figure 6. Maximum-force condition assumed in present model. 


no longer the maximum forces, which occur instead at the vortex locations. For a = 30° shown in fig- 
ure 6(d), for example, the lateral force coefficient at each leading edge of the cross section is equal to only 
0.11, whereas the maximum lateral force coefficient occurs at each primary vortex center and is assumed 
to have a value of 0.39 (the value of (tt/ 2) sin 2 a at a =30°). Setting the total of the lateral force for the 
two primary vortices equal to the leading-edge force for attached flow gives the condition expressed by 
equation (9). 

The steps in applying this boundary condition are included in the bottom part of the flowchart in fig- 
ure 3(a). The vortex center locations for the symmetric case are found by iterating along a curve for a 
given tan a/tan 8, such as one of those shown in figure 5, until an x\,y\ location is found for which equa- 
tion (9) can be satisfied. The resulting curve, labeled “Vortex center locations’’ in figure 7, represents the 
locus of vortex centers for which all three of the previously described boundary conditions are satisfied 
simultaneously. 
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Figure 7. Calculated line of vortex center locations for symmetric vortex flow. 

Determination of Static Vortex Positions and Strengths at Nonzero Roll Angle 

When the wing is rolled, the primary vortices become asymmetrically aligned and have different cir- 
culation strengths (see fig. 8); consequently, all six unknowns in the equation for the complex potential 
(eq. (2)) must be solved for instead of only three. In the present method these unknowns were determined 
by applying the six asymmetric boundary conditions (AS-l-AS-6) stated and discussed in the following 
sections. The first three conditions involve assumptions about how the primary vortices position them- 
selves in response to changes taking place in the underlying potential-flow field when the wing is rolled; 
the remaining three consist of radial-velocity and momentum conditions similar to those applied in the 
symmetric case. 

Asymmetric Boundary Conditions for Wing at Nonzero Roll Angle 

AS-1. Vortex 1 remains on the same attached-flow streamline, deflecting laterally with the streamline, 
as the wing is rolled. That is, 


(26) 

AS-2. Vortex 2 remains on the same attached-flow streamline, deflecting laterally with the streamline, 
as the wing is rolled. That is, 


W 2 =W 2 


<j)= 0° 


(27) 
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Figure 8. Asymmetric vortex arrangement with wing at nonzero roll angle. 

AS-3. At any given roll angle, the two primary vortices both lie on one of the velocity-potential lines 
defined for the body without vortex flow at a zero roll angle. In the a-plane where the cross-flow velocity 
is U* sin a, these lines are defined by (ref. 1 8) 


O ( | )= 0 o = -Uao sinav,. 


l r 2\ 

l + \ 


(28) 


In the t-plane. this condition is met if the two primary vortices lie along the same % = Constant line; that 
is, if 


il=i 2 (29) 

AS-4. The velocity at the center of vortex 1 is in the radial direction. 

AS-5. The velocity at the center of vortex 2 is in the radial direction. 
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AS-6. The force oil the vortex system, in the direction normal to the wing cross section, is equal to the 
rate of change of momentum of the vortex system, in the direction normal to the wing. This condition is 
expressed by the equation: 
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(30) 


Discussion of Asymmetric Boundary Conditions 

Attached-Flow Streamline Conditions (AS-1 and AS-2). These two conditions are based on the 
assumption that the changes in vortex position when the wing is rolled reflect changes in the underlying 
potential-flow field when the wing is rolled. Streamlines of this underlying flow field are shown in fig- 
ure 9 for the X=- and cr-planes (appendix A), the two complex planes for which the vortices can be repre- 
sented mathematically as pure logarithmic vortices. 

The surrounding flow field in these planes can be considered to consist of the two primary vortices 
imposed on a basic attached flow — the flow that would exist if there were no separation and no vortex 
formation. In the a-plane, the flow streamlines are defined, in the general case where the free stream is 
inclined at an angle c|) with respect to the V c -axis, by the equation: 


V = -Uoo sina(;t c sin<|) + y c cos (|)) 


1 - 


2 2 
(a ( , cos<j> - y c sin 4>) _ + (x c sin cj> + y c cos<|>)“ 


(31) 


Figures 9(a) and (b) show representative attached-flow streamlines for the zero roll case; rolling the 
wing in the present model amounted to rotating the free-stream velocity and these streamlines through an 
angle <(j with respect to the body-fixed axes, as shown in figures 9(c) and (d). At a zero roll, the primary 
vortices are located along the 'P; = Constant and hC |V = Constant streamlines, which are shown as 

dashed lines in figure 9. Conditions AS-1 and AS-2 assume that these primary vortices will stay on these 
streamlines and rotate with them as the wing roll angle changes. For example, calculations made in 
this report show that at a = 30°, vortex 1 lies on the streamline = 0.31, and vortex 2 lies on the 

streamline = -0.31. Thus, at a = 30°, the two primary vortices will continue to lie on these 

respective streamlines regardless of the roll angle. 

Constant Velocity-Potential Condition (AS-3). This condition is a result of the flow field being 
irrotational. When the wing roll angle changes and the primary vortices move with the streamlines as the 
streamlines rotate through the angle cjr (as required by conditions AS-1 and AS-2), they must do so in 
such a way that the vortex system as a whole does not rotate with respect to the body. That this require- 
ment can be expressed in tenns of lines of constant velocity potential can be seen by examining the defi- 
nition of these lines in ‘C- and cr-planes used in the present analysis. 


In the ‘C- plane, the velocity potential for <j> = 0° is given by 

£ 

O(K = 0 o = U^a sin a — 
T a 


(32) 
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and in the a-plane, it is given by 
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Therefore, if in the 'C-plane the center of one of the primary vortices is located on the velocity-potential 
line < t ) /c 1 j )=0 o, then according to equations (32) and (33) the nondimensional distance %kla to this vortex 
center can be written as (see also fig. 9(b)) 
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(34) 


Hence, while in the 2-plane, ^/Ja represents the nondimensional distance to a line of constant non- 
dimensionalized velocity potential, that is, the velocity potential divided by aU x sin a (see fig. 9(a)), in 
the a-plane it represents the distance to the centroid of the vortex pair fonned by a primary vortex and its 
image (see fig. 9(b)). (Note that because the two vortices are of equal and opposite strength, this centroid 
lies midway between the primary vortex and its image.) The centroids for the two vortex pairs are desig- 
nated d\ and ch, respectively, in figure 9. 

Consequently, if the two primary vortices always lie on a common <6^ q: line regardless of the wing 
roll angle — that is, if 


then in the 2-plane, 


Vo- = Vo< 

Il = Ii 

a a 


(35) 


(36) 


and there will be no rotation of the two primary vortices (compare fig. 9(a) with fig. 9(b)). In the a-plane, 
there will be no rotation of the line between d\ and di that connects the centroids of the two vortex pairs. 
That is, the relationship 


x r x r 

<j_ _ c 2 


(37) 


will apply, so that the line between d\ and d 2 will always remain parallel to the /Y c -axis (which 
corresponds to the plane of the wing in the A-plane) as the wing roll angle changes (compare fig. 9(b) 
with fig. 9(d)). Applying condition AS-3 will ensure there is no rotation of the complete vortex system in 
either the 2- or a-planes, where the vortices are represented as pure logarithmic vortices. This condition 
should also ensure irrotationality in the A-plane, where because of the conformal transformation involved, 
the vortices can no longer be considered pure logarithmic vortices. 
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Radial-Velocity Conditions (AS-4 and AS-5). AS-4 and AS-5 are the same conditions as condition 
S-l for the symmetric case, but are applied separately at the centers of the two asymmetrically aligned 
primary vortices. Applying these conditions gives two equations similar to equation (13) for the symmet- 
ric case, but each equation contains the circulation strengths of the two vortices. These two equations can 
then be solved simultaneously to determine the two strengths. The procedure for solving these equations 
is described in reference 18. 

Momentum Condition (AS-6). Boundary condition AS-6 is the same as condition S-2 for the sym- 
metric case, except that it is applied here to asymmetrically aligned vortices. Equation (30) represents the 
more general form of equation of motion for the vortex system, containing the coordinates and circulation 
strengths of two vortices instead of just one vortex. 

Application of Nonzero Roll Boundary Conditions 

The basic mechanism by which the two primary vortices move with respect to the wing when the wing 
is rolled is that they move with the underlying potential flow streamlines as these streamlines change 
directions with respect to the wing (conditions AS-1 and AS-2); furthermore, they must both always 
remain along one of the velocity-potential lines that is fixed with respect to the wing (condition AS-3). 
The procedure for determining the static vortex positions and strengths at any given roll angle that result 
from this type of movement is that outlined in the flowchart of figure 3(b). 

The procedure assumes that for a given wing at a given angle of attack and roll angle, the streamlines 
for ()> = 0° have been obtained by applying the symmetric boundary conditions. It begins in the u-plane 
(see fig. 9(a)) where the '%= Constant lines, which correspond to constant <j> = 0° velocity-potential lines 
and do not change with roll angle, are straight lines running perpendicular to the real axis. The stream- 
lines, which do change with roll angle, are defined by 




( X 2 ) 


W = aU x sina 

— cos d) + Im 

-1 

sin 4> 


a 

\\ a 

a 



\ 1 

) 



A | = Constant line, along which the two primary vortices are assumed to lie, is chosen in the "C-plane, 
and by iterating on i] at this value of q using equation (38), a value of 'F can be found that matches the 
H* value for the symmetric case. This process is completed for both primary vortices and establishes 
tentative q and i| coordinates for the primary vortex centers. The vortex circulations can be determined 
by transferring these coordinates to the X-plane and applying the radial-velocity conditions AS-4 and 
AS-5. The momentum condition AS-6 (eq. (30)) is then applied, and if the force on the vortices is not 
equal to the rate of change of momentum of the vortex system, then another value of q is selected and the 
process is repeated. By iterating on q in this way, values of qi,qi and § 2,112 can be found for which 
conditions AS-3 through AS-6 are satisfied, and from these values the vortex center coordinates x\,y\ and 
,Y2,>'2 can be determined in the real X-plane. 

The vortex positions in all three planes considered in the present analysis that result from applying this 
procedure are shown for wing roll angles of 0°, 22.5°, and 45° in figure 10. The x- and y-coordinates of 
the vortex centers determined by the previously described method are given in table 1 for wing roll angles 
from -52.5° to +52.5° and angles of attack from 10° to 40°. 
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(c) (j> = 45°. 

Figure 10. Vortex positions with roll angle shown in three complex planes used in present analysis. 


Determination of Dynamic Vortex Positions and Strengths 

The experimental data indicate that a hysteretic effect on vortex position probably causes wing rock 
(see, e.g., refs. 1-4). When the wing is given a rolling velocity, the primary vortices move away from 
their static positions at a given roll angle; for a positive rolling velocity, the vortex on the side of the 
down-going wing will move away from the wing and above its static position, while the vortex on the 
side of the up-going wing will move toward the wing and below its static position. These displacements 
cause incremental rolling moments, which, at the lower roll angles at least, are destabilizing. In the 
present method, the hysteretic deflections are modeled as displacements of the primary vortices from their 
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static positions; the displacements are determined from the strictly kinematic relationship derived in 
appendix C. Equation (C8) of appendix C gives the displacement in the x-direction as 
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tana 
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, 2U X sin a y 

^2tan8, 


The hysteretic displacement in the y-direction was assumed to be zero, an assumption that is consistent 
with the experimental data of references 1 and 4. 

As derived in appendix C, the hysteretic deflection is equal to the product of the nondimensional lat- 
eral coordinate of the vortex center, the nondimensional roll rate of the wing, and a nondimensional time 
that is constant for the wing and can be expressed in terms of the quantity tan a/ tan 5. Equation (C8) 
assumes each primary vortex is approaching the wing at a velocity of (j) while the wing remains station- 
ary. Using equation (C8) to compute hysteretic effects eliminates the need to employ a source-sink 
distribution on the cylinder surface in the a-plane (the procedure used in refs. 13 and 14 to simulate 
unsteady flow conditions) and greatly simplifies the computational procedure. 

The strengths of the primary vortices in their hysteretic positions were computed by applying the 
radial-velocity conditions at each vortex center as described previously. 

Determination of Wing Rolling Moments 

Calculating the wing rolling moments was a problem in developing the present method because the 
model, which uses just the two primary vortices to represent the vortex flow field, neglects the secondary 
and tertiary vortices that, in a real (viscous) flow, form near the wing leading edges. In a real flow, these 
secondary and tertiary vortices help smooth out the flow at the leading edges so that a Kutta condition can 
be met. Without these vortices, calculated velocities and pressures at the leading edges are infinite, so 
that Blasius integrations around the wing to obtain wing rolling moments do not give the correct values. 
To overcome this problem, methods that did not rely on integrations around the wing had to be used to 
compute the wing rolling moments. 

Static Rolling Moments 

The static rolling moments in this report were computed by using equation (D2) given in appendix D, 
but with the equation modified by multiplying the right-hand side by a constant k. This modified equa- 
tion gives the static rolling moment coefficient for a given wing angle of attack and roll angle as 


C/ = -k-^-sin 2 a sincj) coscj) 


(40) 


In this equation, k is a constant used to account for the difference between the calculated attached-flow 
rolling moments and the experimental moments and is thought to be necessary to account for wing 
thickness effects (as is discussed further in the sections on wing rock results). 
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Moments Caused by Hysteretic Deflections of Vortices 

The wing rolling moment increments resulting from hysteresis in the present method were calculated 
by the method outlined in appendix D, which involved Blasius integrations around the two primary vor- 
tices instead of around the wing. The fonnula used to compute hysteretic moments was 

Cl h =kyAC lh (41) 


In this equation, k is the same constant used in calculating static rolling moments, as described previ- 
ously, and accounts for the difference between calculated and experimental rolling moments. The 
parameter y in equation (41) is a reduction factor used to reduce the rolling moments obtained by inte- 
grating around the vortices to the level of those developed by the wing, as described in appendix D. This 
parameter is a function of tan a/tan 8 (fig. D6) and could be approximated by 

y = 0.00024962 + 0. 10 1 9(tan a /tan 8) - 0.0 1 9932(tan a /tan 8) 2 + 0.00 1 942(tan a /tan 8) 3 (42) 


The incremental coefficient A Ci h in equation (41) is the incremental moment caused by hysteresis in 

the vortex position and was obtained by using the general Blasius equation for moments (ref. 1 9), which, 
in coefficient form, can be written for the wing rolling moment as 


C, = — Re 
6 



A 

a 


' 1 dW 

^f/ooSina df 



(43) 


The incremental moments were obtained by performing two Blasius integrations: one with the primary 
vortices in their static positions x\,y\ at a given wing roll angle, and the second with them in their dis- 
placed positions x\ + A Xk,yk at a given roll angle and roll rate. The incremental moment due to hysteresis 
was then the difference between these two values: 


A C h = Cf 


x k+Axk,yic 



XkJk 


(44) 


Attached-Flow Roll Damping 

Just as a basic attached-flow lift must be added to the vortex lift to get the total lift on a delta wing 
with vortex flow, there is similarly a basic attached-flow roll damping associated with a wing undergoing 
a roll velocity. The theorectial value for the attached-flow roll damping coefficient, as derived by Ribner 
in reference 24, is 


C, 


nA 

32 


(45) 


In the wing rock calculations in the present paper, this attached-flow damping was added to that caused 
by hysteretic effects on vortex position to get the total damping due to roll. 
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Wing Rock Calculations 

To simulate wing rock oscillations, the vortex positions, strengths, and wing rolling moments, 
obtained as described previously, were used in the single-degree-of-freedom equation of motion for the 
delta wing about its roll axis: 

7,^j + l p (j> + /fri ct i on + l h +l s = 0 (46) 


This equation can be written in nondimensional fonn as 


^- + C t -^- + C,. 
qSb p 2U X frlctlon h 


+c h+ c,= 0 


(47) 


In this equation, C/ p is the attached-flow roll damping coefficient, the damping for the wing without 
vortex flow given by equation (45). The symbol C/ ftiction is the rolling moment coefficient due to friction 
in the roll -oscillation apparatus used in a particular wind tunnel test, and C/ A is the increment in rolling 

moment coefficient caused by hysteretic displacements of the vortices, described in the section “Moments 
Caused by Hysteretic Deflections of Vortices.” The symbol Ci s is the static rolling moment, given by the 
suction analogy as 


C[ s = -k —sin 2 a sin <1> cos (]> 


(48) 


To solve this equation, it was put into state variable form by setting 


x \ = <l> 
J 2 =4> 


(49) 


This form allowed the second-order differential equation to be expressed as two first-order equations: 


*1 = x 2 

( 

x 2 = ~ 

V 


c, ^ + c, . . +c, +c, 

l p 2(J friction L h L s 


qSb 


(50) 


In the present method, these equations were subject to the constraints: 


x k + Ax k x k ( y k 

x 2 b 

tana 

ci ci a 

^C/ooSina / 

k 2 tan 5 , 


yk_ = yk_ 

Cl Cl 


(51) 


which are the equations specifying the hysteretic positions of the vortices. 

These equations were solved using a Runge-Kutta scheme, which went through the following steps in 
generating oscillation time histories: 
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1 . 

2. 

3. 

4. 

5. 


The procedure was started at t = 0, cj> = 0, and with an initial value of <j) (1° in most of the calcula- 
tions made in this report) 

(|) was computed from equation (47) 

(|) was integrated over an incremental time At to get (j>, and then (j) was integrated over At to get (j) 
The quantities x k , y k , C k , C/ , and C/l were determined at the new value of <)) 

s ' X k^k 


The new value of (j) was used to compute the product of 

Ax k 


Ci ^ and the hysteretic deflection 
p 2 U x 


6. Values of C k were determined with the vortex centers at x k + A x k ,y k , and those values were used to 

compute C/l . and Ci 

l \x k +Ax k ,y k l h 

7. The time was then set to t + At, and the sequence repeated until a complete time history was 
obtained 


The static vortex positions x k and y k were fed into the equations in the form of cubic-spline curve fits 
to the data of tables 1 (a)-(h). These curve fits gave x k and y k as functions of 4> and were provided for each 
of the angles of attack covered by the tables. The vortex circulation strengths for each primary vortex 
were calculated at each time step by assuming a radial velocity at each vortex center (conditions AS-4 and 
AS-5). The calculations were made for ranges of angle of attack, free-stream velocity, wing sweep, and 
moment of inertia that covered those of the available experimental data. 

Results and Discussion 


Static Vortex Positions and Strengths for Zero Roll Angle 

As shown in figure 11, static vortex positions at zero roll angle are generally in good agreement with 
experimental data and Navier-Stokes calculations (from refs. 2, 1 1, and 25-29): except at the lower val- 
ues of tan a/tan 8, the experimental vortex centers line up along the line of vortex centers determined by 
the present method. The experimental data are for wings with different sweeps tested over different 
ranges of angle of attack. Figure 12 shows that the agreement with experimental data for the present 
method is much better than with either the Brown and Michael (ref. 16) or the Legendre (ref. 30) meth- 
ods, which place the vortex centers much farther outboard than the present method. 

Figure 13 shows the calculated and experimental (from refs. 2, 25, 26, and 27) y- and x-locations as 
functions of tan a/tan 8. At the lower values of tana/tan 8, the vortex center locations determined by 
the present method are farther outboard and away from the wing than the experimental values. These 
differences between calculation and experiment could probably be attributed to the fact that the present 
method neglects the secondary and tertiary vortices which, in the real flow, would fill the space between 
the wing and the primary vortices. At the higher values of tan a/tan 8, where the primary vortices 
become dominant, the agreement between calculation and experiment is quite good. The experimental 
data indicate that wing rock would be expected to occur at values of tan a/tan 8 > 2.0. 
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Figure 1 1. Experimental and calculated vortex center locations. 
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Figure 12. Present method. Brown and Michael, and Legendre methods. 
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(b) x\!a. 

Figure 13. Calculated and experimental x- and y-locations as functions of tan a/tan 8. 


The present calculations showed that the circulation for delta wings could be defined by the single 
curve of Q versus tan a/tan 8 shown in figure 14 (values of Q- horn this curve are presented in table 2). 
Although the vortex circulation T varies directly with distance along a delta wing, for an assumed conical 
flow the value of Q (defined as Ck = 1 72 tm U :tJ sin a) will be the same at each cross section, and therefore 
is a constant for the wing. The curve in figure 14 allows the value of Ck to be determined for any given 
wing at any given angle of attack. 

Figure 14 presents the calculated values of Ck with experimental values obtained from reference 31 
and shows that values from reference 31 are higher over most of the tan a/tan 8 range. However, the 
experimental values in figure 14 were obtained by using a cross-wire anemometer, which is an intrusive 
technique. These values represent the maximum measured in the tests and may contain the effects of 
secondary and tertiary vortices. 
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Figure 14. Calculated and experimental variations of nondimensional vortex strength with tan a/tan S; 4> = 0°. 

As shown in figure 15, the agreement was much better with the data of reference 32, in which circula- 
tion was measured using a novel nonintrusive ultrasound technique. The technique of reference 32 
involves sending ultrasonic pulses in clockwise and counterclockwise directions around a fixed path 
perpendicular to the delta-wing root chord and then determining the circulation from the transit-time 
difference between the two pulses. In figure 15, calculated variations are presented for comparison with 
the data of reference 32 in the form of plots of the nondimensional ratio Tk/UooC against angle of attack, 
for different chordwise stations and two wing sweeps: 60° (figs. 1 5(a)— (d)) and 70° (figs. 1 5(e)— (h)). 
The calculated values in these figures were determined from the calculated curve of Q. versus 
tan a/tan 8 presented in figure 14 by using the relationship 

^ = 2nC k — tan 8 sin a (52) 

I ' qqC c 

The hatched lines in figures 1 5(a) (h) denote the approximate angle of attack at which a vortex burst 
occurred at each chordwise station in the experiments. The data show that at each wing station and for 
both wings, the computed values are in very good agreement with experimental values up to and some- 
times beyond the point of vortex breakdown. 

Static Vortex Positions and Strengths for Nonzero Roll Angle 

Figure 16 shows that, at angle of attack of 30°, calculated vortex center locations for different wing 
roll angles were also in good agreement with the experimental data. Figure 16 presents computed center 
locations along with the wind tunnel results obtained by Jun and Nelson in reference 2, for a wing swept 
80° at roll angles up to 45°. Figure 16(a) shows the center locations relative to the axis system fixed rela- 
tive to the wing; figure 16(b) shows center locations relative to a stationary axis system that rotates with 
the wing. These results are for a wing that is rolled through a positive angle (j), for which the vortex on 
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Figure 15. Calculated values of T/U^c and experimental values obtained by using nonintrusive ultrasound tech- 
nique in reference 32. Hatched lines at each chordwise station indicate approximate angle of attack for vortex 
breakdown at that station; <\> = 0°. 
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(a) Vortex center locations in wing-fixed axes system. (b) Wing and vortex center movement relative to fixed 

axes. 

Figure 16. Experimental and calculated vortex center locations for different wing roll angles at a = 30°. 

the right wing panel is located inboard and close to the wing, while that on the left panel is located farther 
outboard and away from the wing. Figure 16 shows the agreement to be very good at roll angles of about 
30° or below. Above about 30°, the calculated results put the right vortex farther inboard and the left 
vortex farther outboard than the experimental results indicate. 

The agreement is not as good at the higher angles of attack. The effect of angle of attack on the com- 
puted positions is shown in figure 1 7, which shows that the center locations move away from the wing as 
the angle of attack increases, but the increase in distance from the wing becomes small at the higher val- 
ues of a. The wind tunnel results, on the other hand, show vortex centers that are farther from the wing at 
a = 45° than the calculated results indicate they should have been. This difference could possibly be 
caused by a bending away from the wing due to the influence of the flow aft of the wing. At a = 40°, the 
wind tunnel results show vortex breakdown, indicating that some of this bending has occurred. 

The calculated effect of positive roll angle on nondimensional vortex strength is shown in figure 18. 
Except for the lowest angle of attack, the circulation is higher for vortex 1 , the vortex closer to the wing. 
This vortex tends to maintain about the same circulation throughout the range of 4>, whereas vortex 2, 
which moves out and away from the wing at positive roll angles, loses circulation with increasing values 
of c|). This type of variation is consistent with the increase in the static restoring moment on the wing that 
occurs when rjj increases and no hysteretic effects are present. 
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Figure 17. Static vortex center locations for various wing roll angles and angles of attack. 

Wing Rock Time Histories 

Comparison of Calculated and Experimental Time Histories 

The results of the wing rock calculations show that, while in general the present method was able to 
model the basic wing rock motion for delta wings, the motions were highly dependent on the static rolling 
moment characteristics of the wing. As will be discussed, when the variation of static C/ versus <]> 
obtained by using the leading-edge suction analogy was used, the rates of buildup were more rapid and 
the maximum oscillation amplitude and frequencies were higher than in the experiment; however, when 
this variation of static C/ versus (|) was adjusted to bring it into closer agreement with the wind tunnel 
variation, the oscillations were brought into close agreement with the experiment. 

Figure 19 presents calculated and experimental wing rock oscillations for a delta wing swept 80° at an 
angle of attack of 30°. The experimental data are those obtained in reference 4 in free-to-roll tests of a 
wing that was mounted on an air-bearing support system that had essentially zero friction. This wing has a 
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<K deg 

Figure 18. Effect of wing roll angle on nondimensional vortex strength (radial velocity assumed at each primary 
vortex center). 

root-chord length of 1.3976 ft, a span of 0.4921 ft, an area of 0.3439 ft 2 , and a roll axis moment of inertia 
of 6.39 x 10-4 slug-ft 2 . Figure 19 shows both the roll angle time histories and the variations of the static 
rolling moment coefficient with roll angle. The variation of C/ versus tj> for the experimental data was 
determined from the roll angle time history as described in reference 4. 

Figure 19 shows a close tie-in between the roll angle time histories and the static C/ versus (j) curves. 
Both the experimental and the calculated time histories show that frequency generally decreased as 
amplitude of the oscillation increased; this inverse relationship is a consequence of the fact that the oscil- 
lation frequency depends primarily on the slope of the C/ versus (j) curve, and this slope decreases with 
increasing roll angle. 

The calculated time history (fig. 1 9(b)) shows a more rapid buildup to a higher maximum amplitude 
than for the experimental time history (fig. 19(a)) (the calculated time history reaches (|) niax of about 52.2° 
instead of about 49.4° for the experiment) and has a higher frequency at the maximum amplitude (1.43 Hz 
compared with about 1.17 Hz for the experiment). These differences could be attributed to the substantial 
differences in the C/ at a given f|) between calculation and experiment. When the calculated C/ versus (j) 
curve was generated, it was assumed that the static wing rolling moments were those given by equa- 
tion (40), and that the wing had a flat-plate cross section for which k = 1 .0. This assumption gave much 
higher calculated static rolling moments at a given (j> than those shown for the experiment in figure 19(a). 
Since the hysteretic moments were assumed to be proportional to the static moments, the hysteretic loops 
of the calculated rolling moment variations in figure 19(b) were substantially wider that those for the 
experiment in figure 19(a). The higher rolling moment increments produced by these hysteretic loops 
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(a) Reference 4 experiment. 



t, sec Q 


(b) Calculated with assumed static rolling moment reduction factor k = 1.0. 



(c) Calculated with assumed static rolling moment reduction factor k = 0.649. 

Figure 19. Measured and predicted roll angle time histories and Cj variation with <j) at a = 30°, <|) = 1°, and A = 80°. 


were the reason for the faster rate of buildup and higher maximum amplitude of the calculation compared 
with the experiment. The higher oscillation frequencies for the calculation compared with the experiment 
could be attributed to the slopes of the C/ versus (|) curve being higher than those for the experiment at all 
roll angles. 

Like figure 19(b), figure 19(c) shows a calculated oscillation time history but with an assumed value 
of k = 0.649. Figure 19(c) shows that with this new value of k, the calculated time history is brought into 
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close agreement with the experimental data of reference 4 (shown in fig. 19(a)): the calculated rate of 
buildup approximately matches that of the experiment, the calculated maximum oscillation amplitude is 
about 47.5° versus 49.4° for the experiment, and the calculated frequency at maximum amplitude is about 
1.19 Hz versus 1.17 Hz for the experiment. This closer agreement would then indicate that if the static 
rolling moment characteristics of a wing (or possibly a wing-body configuration) could be determined 
fairly accurately, the wing rock motions could be approximated fairly accurately. 

The lower value of C/ at a given (j) for the experiment could possibly be due to a wing thickness effect. 
The calculations assume a thin flat-plate wing; the experimental wing had a thickness of 0.25 in., giving it 
a thickness-to-span ratio at the wing trailing edge of 4.23 percent. This thickness-to-span ratio increases 
as the distance from the wing apex to a given wing cross section decreases and theoretically becomes 
infinite at the apex. The type of beveling at the wing leading edge could also account for differences in 
calculated versus experimental results. Reference 9 shows significant effects of different types of bevel- 
ing on wing rock amplitudes. For example, the wing tested in reference 1 had a sharp upper surface with 
beveling on the bottom surface; according to reference 9 this type of beveling would lower the severity of 
the wing rock. 

Figure 20 provides another set of wing rock time history data for comparison: data from reference 7, 
for a wing swept 80° at an angle of attack of 25°. The wing tested in reference 7 had a root chord of 
1.0827 ft, a span of 0.3806 ft, an area of 0.3054 ft 2 , and a roll axis moment of inertia of 1.566 x 10 -4 
slug-ft 2 , and was tested at a free-stream velocity of 49.22 ft/sec. The oscillation shown was initiated at a 
wing roll angle of 10°. The static rolling moment variation with (j) for this wing was not provided in ref- 
erence 7, nor is it known how much friction was present in the mounting apparatus used in the tests. This 
wing had a thickness-to-span ratio at the trailing edge of 1.72 percent and had its sharp edge at the lower 
surface with beveling on the upper surface. By setting k = 0.87 and by assuming a value for the structural 
present in the tunnel roll apparatus for the tests, the calculated time history could be brought into close 
agreement with the experiment for the wing of reference 7. A value of C/ fti ti = 0.03 was assumed, 
which is about half the value of C/ for the wing with attached flow. 

Comparison of Calculated and Experimental Hysteretic Deflections 

When the approximation to the experimental static rolling moments was used, the calculated hysteretic 
deflections were in good agreement with experimental data of reference 2. Calculated hysteretic deflec- 
tion envelopes that assume k = 0.649 are presented for comparison with experimental results in figure 21, 
which shows the x versus y position of the vortices during a steady state oscillation at maximum ampli- 
tude. Figure 21 shows the calculated envelopes are slightly longer and narrower than the 
experimental envelopes, but overall the agreement is good. The y (lateral) and x (normal) hysteretic 
deflections are presented as functions of (|) in figures 22(a) and 22(b), respectively. In figure 22, the 
experimental envelopes are not quite symmetrical; nevertheless, overall there is good agreement between 
calculation and experiment. 

Causes of Wing Rock and Its Limited Amplitude 

The calculated results of the present method verify that — as has been indicated by the experimental 
results of references 1 through 4 — hysteretic deflections of the leading-edge vortices cause wing rock. 
The results show that the hysteretic deflections at the higher angles of attack and the lower roll angles 
produce the unstable rolling moments that drive the wing oscillation to higher amplitude. Moreover, the 
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t, sec 


(a) Reference 7 experiment. 



t, sec 

(b) Calculated with assumed k= 0.87, Cl. . =0.03. 

v 7 1 friction 

Figure 20. Calculated amplitude time history and experimental data of reference 7 at a = 25°. 

calculated results indicate the source of the restoring moments that limit the amplitudes of the oscillation. 
When the wing is rolled, the asymmetric vortex arrangement consists of one vortex that is located farther 
outboard and away from the wing than the other vortex, which is inboard and close to the wing. The 
restoring moments can be attributed to the movement of the outboard vortex to points farther away from 
the wing and an increase in the strength of the inboard vortex. 

The rolling moments caused by wing rolling velocity are shown for angles of attack of 10°, 20°, and 
30° in figures 23(a), (b), and (c), respectively. Each figure shows total roll damping as a function of non- 
dimensional roll rate for wing roll angles of 0°, 15°, 30°, and 45°. The total roll damping shown is the 
sum of the attached-flow roll damping (equal to the attached-flow roll damping coefficient times the non- 
dimensional roll rate) and A C/ , the incremental hysteretic damping. 
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Present method 



Figure 21. Calculated and experimental normal hysteretic displacement of primary vortices as function of lateral 
displacement at a = 30°. 


Figure 23(a) shows that at a = 10°, the roll damping curves for each wing roll angle fell along a single 
straight line; this colinearity was the result of hysteretic vortex deflections at this angle of attack being 
very small, so that the hysteretic roll damping moments were essentially zero. The slope of these curves 
was approximately equal to C/ — the roll damping coefficient for the wing with attached flow — with a 

negative slope (negative C/ ) indicating stable damping (that is, a positive roll rate would produce a 
negative moment that opposes the wing motion). Because the roll damping at a= 10° was stable, wing 
rock did not occur at this angle of attack. 

At a= 20° (fig. 23(b)), hysteretic vortex deflections were larger, and these larger deflections 
decreased the stability to very low levels at c|) = 0° and 15°; however, the damping moments remained 
either neutrally or slightly stable so that, although the wing was on the verge of wing rock at this angle of 
attack, wing rock still did not occur. 

At a = 30° (fig. 23(c)), the rolling moments show a high degree of instability at roll angles of 0° and 
15°. Wing rock did occur at this angle of attack. 

At all three angles of attack, the roll damping moment variation was stable at a roll angle of 30°, 
where it was approximately equal to the attached-flow roll damping coefficient; at <x = 20° and 30°, it 
became more stable when the roll angle was increased to 45°. The sketches in figures 23(b) and (c) show 
that at least part of this stable damping at the higher roll angles could be related to the movement of the 
vortex on the upward-deflected wing panel to points outboard and farther away from the wingtip. At 
values of c|) less than 30° both vortices lie inboard of the wingtip, but at (j) = 30° the vortex on the upward- 
deflected panel is located slightly outboard of the wingtip and continues to move farther outboard and 
away from the tip as the roll angle increases. Much of the stable damping could be attributed to this vor- 
tex moving outboard of the wingtip to where it had much less of an effect on the wing. 

At the same time, the calculated circulations, such as those presented in figure 24, show that the 
inboard vortex in these sketches would have a higher circulation (C 2 in fig. 24 would be more negative) 
when the wing roll angle is positive and increasing. Although the inboard vortex is farther from the wing 
than it would be in its static position, it produces a positive rolling moment that acts to restore the roll 
angle to zero. These two effects combined — the lessening effect of the outboard vortex and the higher 
circulation of the inboard vortex — produce the positive restoring moments that limit the amplitudes of 
wing rock oscillations. 
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Figure 24. Variation of nondimensional vortex strengths of primary vortices with wing roll angle for wing rock 
oscillation at a = 30°. 


Effect of System Parameters on Oscillation Amplitude and Frequency 

The calculated results showed that for wings with a given sweep and cross section shape, the wing 
rock oscillations could be defined in tenns of two parameters that were functions of wing angle of attack 
only. These parameters were the maximum amplitude of the oscillation <j> max and the reduced fre- 
quency f-yjli / qSb . These parameters were affected by changes in wing sweep and changes in k, the wing 
cross section shape parameter (see eq. (40)); also, in wind tunnel tests they were affected by the amount 
of friction in the mounting apparatus. But for wings with the same sweep, cross section, and test condi- 
tions, the variation with angle of attack would be the same, regardless of the free-stream velocity, the 
wing size, or the wing moment of inertia. 

Calculated variations of (j> max with a are presented for comparison with experimental variations for a 
delta wing with 80° and 75° sweep in figures 25(a) and (b), respectively. Figure 25(a) shows wing rock 
oscillations for experimental data that represent wings with different thickness and different types of 
leading-edge beveling. The experimental data of reference 9 show wing rock starting at an angle of attack 
of about 17°, compared with about 22° for the experimental data of reference 4. An abrupt drop-off in 
cjimax at about a = 30° can be seen in the data of reference 4, which, according to reference 1, is due to 
vortex breakdown. Two calculated variations are also shown: one for the wing of reference 4, for which 
it was assumed that k = 0.649 and that no structural damping was present, and the other for the wing of 
reference 7, for which it was assumed that k = 0.87 and that there was structural damping in the amount 
C/^ction = 0 03 present in the tunnel roll apparatus. Overall, the calculated results agree with the experi- 
mental data. The calculated results with k = 0.649 generally fall between the experimental data of refer- 
ences 4 and 9, which had similar types of leading-edge beveling; the calculated results with k = 0.87 
and C/ ftiction = 0.03 for the wing of reference 7 are in good agreement with the experimental data of 

reference 7. 
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(a) A = 80°. 



a, deg 


(b) A = 75°. 

Figure 25. Calculated and experimental variations in maximum amplitude of wing rock oscillations with angle of 
attack. 


46 



By comparing both calculated and experimental data from figure 25(a) (A = 80°) with that of fig- 
ure 25(b) (A = 75°), one can see that reducing the wing sweep reduced the maximum amplitude of a wing 
rock oscillation at a given angle of attack. This reduction in (j) m ax is the result of the fact that, although 
wings with different sweeps have the same static C/ at a given angle of attack and roll angle (since they 
were computed by using eq. (40)), the wing with the lower sweep has a lower value of tan a/tan 8 at a 
given angle of attack. With a lower value of tan a/tan 8, the vortex circulations and hysteretic displace- 
ments are lower at a given roll angle; hence, the chiving moments are smaller and the wing rock oscilla- 
tion does not reach as high a maximum amplitude. 

The effects of angle of attack on wing rock reduced frequencies are shown in figures 26 and 27, where 
the calculated results from the present method are presented along with experimental results from refer- 
ence 7. The data from reference 7 are the only set available for demonstrating the effects of free-stream 
velocity, wing sweep, and wing inertia on wing rock characteristics. In reference 7 reduced frequency is 
presented in the form fb/U x , which makes it independent of free-stream velocity. 

Figure 26 presents calculated and measured values of fb!U x for two wing sweeps and different wing 
moments of inertia. Calculated values of fb/Uao increased with angle of attack for both wing sweeps. For 
the wing swept 80°, the experimental values of the reduced frequency remained approximately constant 
with angle of attack below a =40° for the low wing inertias. At angles of attack above 40°, the experi- 
mental values tended to increase and come into closer agreement with the calculated values as the wing 
inertia increased. For the wing swept 75°, the experimental values of reduced frequency increased with a 
and were in good agreement with the calculated values for both wing inertias shown. 

As shown in figure 27, the calculated curves for different wing inertias from figure 26 could be 
reduced to just two curves Tor wing sweeps of 80° and 75° — by expressing reduced frequency in 
the form f^TjJqSb. The reason for this is that, for a spring- and-mass system such as that defined by the 
equation of motion for the wing rock model (eq. (47)), the oscillation frequency is primarily determined 
by the spring-constant term, that is 


/- 


h 


(53) 


Since, for wings with the same cross section shape, the Cy versus ()> variation depends only on angle 
of attack (see eq. (40)), the reduced frequency fyjI[/qSb should depend only on angle of attack and be 
independent of wing size and inertia and free-stream dynamic pressure. 

The experimental values of f^lJJqSb plotted in figure 27 show that the data for the two lower iner- 
tias for the wing swept 80° did not collapse into a single curve, and the values were lower than the calcu- 
lated values; however, as with the velocity-reduced frequency fblU<x (shown in fig. 26), the data did tend 
to collapse and the reduced frequencies came into closer agreement with the present calculations as the 
inertia was increased. For the higher wing inertias for the wing swept 80°, and for all inertias for 
the wing swept 75°, the calculated reduced frequencies of the form f\jl / / qSb were in generally good 

agreement with the experimental values. 
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(a) A =80°. 



Figure 26. Calculated and experimental variations with angle of attack of reduced frequency fbUJaa. 
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Exp. 
(ref. 7) 



(a) A = 80°. 



(b) A = 75°. 

Figure 27. Calculated and experimental variations with angle of attack of reduced frequency f-JTJJqSb . 


Conclusions 

A discrete-vortex method has been developed to model the vortex flow over a delta wing and investi- 
gate the wing rock problem associated with these wings. The results indicate the following conclusions: 

1. Static vortex positions at zero roll angle were in good agreement with experimental results. The 
present method gave much better agreement with experimental data than either the Brown and 
Michael or Legendre methods, which placed the vortex centers much farther outboard than the 
present method. 
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2. For angles of attack below those for which vortex breakdown occurred, calculated vortex circula- 
tions were in excellent agreement with those measured by Johari and Moreira in “Direct Measure- 
ment of Delta-Wing Vortex Circulation” (AIAA Journal, December 1998, pp. 2195-2203) using a 
nonintrusive ultrasound technique. 

3. Static vortex center locations at nonzero roll angles were also in good agreement with experimental 
data. At angles of attack a of 35° and 45°, the wind tunnel results showed vortex centers that were 
farther from the wing than those of the calculated results. This difference may have been caused 
by the bending away from the wing due to the influence of the flow aft of the wing. At both 
a = 35° and 40°, the wind tunnel results showed vortex breakdown, indicating some of this bend- 
ing had occurred. 

4. The present method was able to model the basic features of wing rock. The rates of buildup of the 
oscillations were generally more rapid and the maximum oscillation amplitude and frequencies 
were higher than the experiment when the variation of the static rolling moment Cj versus roll 
angle <j) obtained from the leading-edge suction analogy was used. The oscillations were brought 
into very close agreement, however, when an approximation to the variation of the static C/ versus 
cj) obtained from the wind tunnel results was used. 

5. When the approximate variation of static C/ versus (|) obtained from wind tunnel results was used in 
the calculations, the calculated hysteretic deflections were overall in good agreement with the 
experiment, indicating that a strictly kinematic relationship — derived by using the time analogy 
and conical-flow assumptions — could be used to calculate these deflections. 

6. The present method was able to verify the cause of wing rock and indicated the reason for the 
limited amplitude. At roll angles of about 30°, a lag in vortex position produced a destabilizing 
rolling moment increment; above about 30°, it produced a stabilizing increment. This shift from 
destabilizing to stabilizing rolling moments was the result of one vortex moving outboard of the 
wingtip and decreasing in strength at wing roll angles above about 30°. 

7. The calculated results showed that for wings with a given sweep and cross section shape, the 
maximum amplitude and the reduced frequency fyfTjJqSb would be functions of wing angle of 
attack only. Calculated variations with angle of attack of the maximum amplitude of the wing rock 
oscillations generally agreed with the experimental data. Both the experimental and calculated 
results showed that reducing the wing sweep reduced the maximum amplitude of the wing rock 
oscillation for a given angle of attack. 

8. For wings swept 80°, experimental values of f\jl / / qSb were generally lower than the calculated 
values for the lower wing inertias, but came into closer agreement with the calculations as the 
inertia was increased. For the higher wing inertias, calculated and experimental reduced frequen- 
cies were in good agreement for wings swept both 80° and 75°. 
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Appendix A 

Transformations for Complex Planes Used in Present Analysis 


The transformations equations presented in this appendix allow transformation back and forth between 
any of the three complex planes, shown in figure A 1 , that were used in developing the present method. 

In the (1-plane (fig. Al(a)), the flow at zero roll angle is parallel to the flat plate representing a delta 
wing cross section. This plane is used in the Brown and Michael analysis of reference 16, wherein it 
serves as the primary complex modeling plane in which the problem was set up before transforming to the 
physical plane. It could not be used in this capacity in developing the present method, however, because 
it did not give the correct velocities at the vortex centers in the physical plane, as required by the radial- 
velocity boundary conditions of the present model. To determine the velocity at a given vortex center, the 
effects of that vortex on itself must be neglected. But in the 'C-planc one vortex also serves as the image 
of the other, so that two singularities have to be dealt with at a given center. Setting the velocities at both 
singularities to zero gives incorrect values for velocity at the vortex center in the physical plane. 

In the present method the 'C-plane is used primarily to define the vortex locations when the wing is at a 
nonzero roll angle. In this plane, lines of constant velocity potential at a zero roll angle correspond to 
| = Constant lines and so are easily defined. The procedure used to iterate on | in order to determine 
vortex positions for nonzero roll angles is described in the section “Application of Nonzero Roll Bound- 
ary Conditions.” 

The a-plane (fig. A 1(b)) is the basic modeling plane used to generate flow fields in the t- and 
A-planes. In this plane, the vortex flow is represented by pairs of simple logarithmic vortices (see eq. (2)) 
for which the momentum is known. The images placed inside the cylinder ensure tangent flow at the 
surface of the body in all three planes. Since the images are separate from the primary vortices, the 
velocities calculated at the primary vortex centers (using eq. (4)) will have the correct values, so that the 
radial-flow velocity condition can be applied. 

The A-plane (fig. Al(c)) is the physical plane in which the actual flow about the wing cross section is 
represented, and is the one in which all boundary conditions must be satisfied. 


Points in one of these complex planes can be transformed into points in any of the others by using one 
of the following transformations: 


t = -yjx 2 + 1 

Y ° 2 

C = o + — 
4 o 


A = o 

4 o 


r-Vri-i 


a 


a = 


A ±y/A. 2 +a 2 


(Al) 
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In nondimensional form, these become 
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Appendix B 

Description of Forces Acting on Delta Wing 

The description of delta wing forces in this appendix is intended to clear up some of the confusion 
existing in the literature about the nature of these forces, particularly about the force on the feeding sheet 
considered in previous discrete-vortex models such as those of references 13, 14, 16, and 17. 

According to the present model, three forces act on the delta wing: one steady and two 
unsteady. These forces can be computed by performing Blasius integrations (see refs. 19 and 20) around 
a closed path that encloses only the wing in the physical /.-plane. Steady forces are computed using the 
formula: 

(b » 

Unsteady forces are computed by using the time-dependent formula: 

F x -iF Y =-i P— (6 WdX (B2) 

■^unsteady -'unsteady dt 

In these equations, W is the complex potential for the flow (see eq. (2)), and can be written in terms of 
coordinates in the /.-plane as 



The unsteady forces arise from the fact that this complex potential is a function of time because the 
wing semispan a and the vortex circulation I)- vary with time. That is, both quantities vary with distance 
along the three-dimensional delta wing and, therefore, according to the time analogy, vary with time in an 
equivalent two-dimensional time-dependent plane. The vortex center coordinate 7,/- also varies with dis- 
tance or time, but appears in the complex potential only in its nondimensional form X^/a, which is inde- 
pendent of time. That is, for a given wing at a given angle of attack and roll angle, X^/a is the same for 
each primary vortex at each wing cross section and is therefore independent of both distance along the 
wing and of time. 

The three forces that can be obtained using the steady and unsteady Blasius formulas are discussed in 
the following sections. 

Static Pressure Force 

The static pressure force is obtained by considering all quantities in equation (B3) to be constant and 
independent of time, and therefore it is considered a steady force that can be determined by using equa- 
tion (Bl). In deriving this force, the primary vortices are assumed to have fixed positions and to have 
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fixed circulations; this force then represents the static pressure force that the vortex system exerts on the 
wing as a result of its inducing high velocities and therefore high pressures on the wing. This force can 
be written in coefficient form as 


C 


x ' 


iCy = — CD 


c. 
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dW 


C/ooSina dk 


\ 2 

(k\ 

d \ 

t) 


(B4) 


The opposite of this force is the force that the wing exerts on the vortex system, and it is this static 
pressure force on the vortex system that causes the primary vortices to change their positions and hence 
their momentum, as expressed by the boundary conditions S-2 and AS-6. 


Attached-Flow Force 


The attached-flow force is obtained by holding 1’/. and k^ constant when taking the time derivative in 
equation (B3) while allowing the local wing semispan a to vary with time. This is the basic force that the 
wing would develop if the flow were attached and there were no vortex formation. Although the 
attached-flow force is usually considered to be a steady force, there would be no such force in potential 
flow if the wing cross section were not expanding. In the context of the present model, it is the force that 
results strictly from the fact that the wing cross section expands with distance along the delta wing or with 
time in an equivalent two-dimensional flow. 

The integral involved in obtaining the attached-flow force can be written in terns of the wing semi- 
apex angle and angle of attack using the time analogy. The derivative with respect to time in equa- 
tion (B2) can be written as 


d _ d da dl 
dt da dl dt 


— Una cos a — = tan 6 U x cos a — 
/ da da 


(B5) 


With the time derivative written this way, the attached-flow force, in coefficient fonn, becomes 
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(B6) 


where the integration is to be carried out around any path enclosing the wing. 


It can be shown that the lift coefficient computed for the wing by using equation (B6) is 
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(B7) 


which, when the small-angle approximation for sin a is made, becomes the formula derived by Jones for 
slender delta wings in reference 33. 

Force Caused by Changes in Vortex Circulation 

This force is obtained by holding the wing semispan a constant and allowing T to vary in equa- 
tion (B2). It is the force that arises from the increases in T that take place as the primary vortices grow 
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with time in the two-dimensional time-dependent plane, or with distance along the three-dimensional 
delta wing. This force is the same as that on the straight-line feeding sheet that serves as a branch cut in 
the Brown and Michael (ref. 16) and similar models. That is, these previous models assume a force on 
the branch cut of 


(B8) 

at 

This force then is equal to the unsteady pressure pdT/dt times the length of the branch cut. 

It should be noted, however, that the branch cut for which this force is calculated cannot be taken 
arbitrarily; more specifically, it cannot be taken to lie between the primary vortex and a separation point 
at the wingtip as assumed in the Brown and Michael model. The logarithmic functions in equation (B3) 
have branch points z 0 and at the image point r~ la 0 , thus the branch cut needed to make them single - 
valued should lie between these two points. In the a-plane, this cut is along the straight radial line 
between the primary vortex and its image; in the /.-plane, it becomes a curved line extending from the 
wing to the vortex center, as shown in figure 2(b). To calculate the force using equation (B3), the inte- 
gration must be carried out around a cut that connects the vortex and image branch points. An integration 
carried out around a path enclosing only the feeding sheet of the Brown and Michael model will give zero 
force, rather than the force given by equation (B8). 

Although the force given by equation (B8) is to be interpreted mathematically as a pressure acting 
over a certain length in the flow, it can be given no such physical interpretation in the real flow. It is 
strictly the unsteady pressure force that results from changes in the strength of the vortex system. There 
is no physically visible feeding sheet in the real flow along which unsteady pressures act. 

As discussed in the main text in the section “Discussion of the Symmetric Boundary Conditions: 
X-Axis Momentum Condition,” this force is neglected in considering the force-momentum balance on the 
vortex system because the unsteady force caused by changes in vortex circulation strength is exactly 
equal to the rate of change of momentum caused by changes in vortex circulation, regardless of the loca- 
tion of the vortex center. No new information is obtained by equating this unsteady force to the rate of 
change of momentum caused by changes in circulation; therefore, this force does not play a roll in deter- 
mining vortex position and strength. 

Total Delta Wing Force 

Because secondary and tertiary vortices are not included in the present model, adding the steady and 
unsteady forces described previously will not give the correct value for the total force on the wing. The 
total normal force for a delta wing, however, can be closely approximated by using the Polhamus leading- 
edge-suction analogy (ref. 23), which gives a normal force coefficient due to the vortices (for the general 
case with the wing at a nonzero roll angle) of 

C)y = jt sin 2 a cos 2 (j) (B9) 
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Appendix C 

Derivation of Equation for Vortex Hysteretic Deflections Caused by Wing 
Rolling Motion 

The vortices in the wake of a delta wing do not respond instantaneously to wing motion and therefore 
are subject to hysteretic deflections as the wing rotates about its longitudinal axis. The available experi- 
mental data show very little, if any, hysteretic deflection in the lateral direction (the data of ref. 4, for 
example, show v-positions of the vortices to be approximately the same as the static positions at any time 
during a rolling oscillation). The deflections in the x-direction (normal to the wing), however, can be 
substantial. The model in the present method makes use of the two-dimensional time analogy and also 
assumes that the primary vortices deflect in a conical manner as the wing rotates, as shown in figure Cl. 
In figure Cl, the line of vortex centers is assumed to deflect as a straight line, with the deflection Ax at 
any station being proportional to the distance from the wing apex, as required if the flow is to remain 
conical. 



According to the two-dimensional time analogy, the flow over a delta wing is analogous to that for a 
two-dimensional plate whose span and the flow around it are changing with time. In the two-dimensional 
flow field, the vortices are increasing in strength and moving radially outward from the center line of the 
plate with time. If the plate is rotating, each vortex will also be moving either toward it or away from it at 
a velocity of y^( j> and will deflect toward or away from the plate through a distance of 

= fotfdt (Cl) 

Multiplying and dividing the right-hand side of equation (Cl) by aU x sin a gives 

Ax>. = — C/qo since f a dr (C2) 

a t/ooSina •'O 

If the flow is conical, the ratio yija in equation (C2) will be a constant, since the vortex positions are 
changing in proportion to changes in the semispan; the variable a by itself, of course, will vary linearly 
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with time. According to the time analogy, time is related to chordwise distance along the three- 
dimensional wing by 


dt = — (C3) 

f/ooCosa 

For the delta wing, local semispan is related to chordwise distance by 

a = 1 tan 8 (C4) 

Making these substitutions into equation (C2) and carrying out the integration over the length of the root 
chord of the three-dimensional wing gives the hysteretic deflection at the wing trailing edge: 

2 

(Avr.) tP = — tana f /tan8r/Z = — tana— tan 8 (C5) 

a f/ooSina a Lisina 2 


For the delta wing, L = b ! 2 tan 8, so that equation (Cl) can be written as 
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The deflection at the base, nondimensionalized in tenns of the wing semispan at the base, then becomes 
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Because the flow is assumed to be conical, both the hysteretic deflection Ax and the wing semispan a vary 
linearly with chordwise distance /. With this type of variation, the ratio of hysteretic deflection to wing 
semispan will be constant along the wing, so that at any cross section this ratio is equal to the ratio at the 
trailing edge. Equation (C7), therefore, gives the nondimensional hysteretic deflection at any cross 
section as 


Ax k (Ax i .) te y k 
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(C8) 


Equation (C8) shows the hysteretic deflection to be equal to the product of ( 1 ) the nondimensional 
lateral coordinate of the vortex center, (2) the nondimensional rolling velocity of the wing, and (3) a 
nondimensional time, which is constant for the wing and can be expressed in terms of the quantity 
tan a/tan 8. 
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Appendix D 

Wing Rolling Moment Calculations 


As stated in the main text, calculating the wing rolling moments was a problem in developing the 
present method because of the simplifying assumptions made. The model neglects the secondary and 
tertiary vortices and does not apply a Kutta condition at the wing leading edges; therefore, the velocities 
and pressures at the leading edges are infinite. Under these conditions, Blasius integrations around the 
wing do not give the correct moments, presenting a problem in calculating both the static wing rolling 
moments and the rolling moments produced on the wing by hysteretic deflections of the primary vortices. 
For the static rolling moments, this problem was solved by assuming that these moments could be closely 
approximated as being those for the wing with attached flow. For the hysteretic moments, it was resolved 
by integrating around the primary vortices instead of the wing to obtain the moments, and then reducing 
these by a constant factor to obtain the wing moments. The equations for the static and hysteretic 
moments are discussed in the following sections. 

Calculation of Static Wing Rolling Moments 

These static moments were determined by integrating around the wing cross section with an assumed 
attached flow. This integration was done by using the Blasius fonnula for moments (see ref. 19 or 20), 
which in coefficient fonn can be written as 


C, = -- Re 
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The complex potential W(k) used in this equation was that for the free-stream flow about the wing with 
no vortex flow (that is, eq. (2) with C\ and C 2 set to zero). The integration was carried out around a path 
that enclosed the complete wing cross section. 

The moment coefficients obtained from equation (Dl) can be written in closed fonn for a wing at roll 
and/or sideslip as follows: 

For a rolled wing with zero sideslip, 


~ 71 • 2 • , , 
C/=-— sin a sin <p cos <j) 


For a sideslipped wing with zero roll, 


C; = sin a sin |3 cos (3 

For a wing that is both rolled and sideslipped, 

2 2 2 1 2 
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Values of rolling moment coefficients calculated by using these relationships are presented with 
experimental results in figures D1 through D3. Figures D1 and D2 show the variations with angle of 
attack for different roll and sideslip angles, respectively; figure D3 shows the variations with sideslip and 
roll angle at angles of attack of 10° and 20°. These figures show that the analytical values are in generally 
good agreement with experimental data for angles of attack below the stall. 

Calculation of Hysteretic Rolling Moments 


Because the attached-flow rolling moments used to approximate the static wing rolling moments in 
equations (D2)-(D4) were functions of wing angle of attack, roll angle, and sideslip angle alone — and not 
vortex position — they could not be used to determine the hysteretic rolling moments when the wing was 
oscillating. These hysteretic moments had to be determined by Blasius integrations made with the 
primary vortices included in the flow field. 

With vortices included, it was found that the rolling moments calculated by Blasius integrations 
around only the wing would not limit the amplitude of the wing rock oscillations. This failure to limit the 
amplitude was believed to be the result of the infinite pressures acting at the wing leading edges in the 
present model. These infinite pressures give a distorted pressure distribution on the wing so that moments 
cannot be computed accurately. 

To work around this problem, the rolling moments used in the calculations were those obtained by 
integrating around only the primary vortices. The moments obtained this way are the result of finite 
forces that act at the centers of the primary vortices to produce a moment about the wing moment center. 
It was believed, therefore, that incremental changes in the magnitudes and positions of these forces 
caused by hysteresis would produce incremental moments that were an accurate indication of those 
affecting the wing. 


As indicated previously, the incremental moments due to hysteresis were obtained by using the 
Blasius equation (Dl), and for W(K), using the full complex potential for the flow with vortices (eq. (2) of 
the text with Ci and C 2 not equal to zero). The integrations were carried out around a path that included 
only the two primary vortices. To determine the hysteretic moments, two Blasius integrations were made: 
one with the primary vortices in their static positions Xk,yk at a given wing roll angle; the second with 
them in their displaced positions Xk + &Xk,yk at a given roll angle and roll rate. The incremental moment 
due to hysteresis was then the difference between these two values: 


c l h = c l 


— c 

A'/.+/Yv£ ,V£ f X k. k. 


(D5) 


As shown in figure D4, the rolling moments obtained by integrating around the primary vortices were 
substantially more negative than those obtained by integrating around the wing only and, therefore, had to 
be reduced to obtain the moment increment that would actually be felt by the wing. This reduction was 
accomplished by multiplying the computed moment increment by reduction factor 7 , which was obtained 
from plots of wing-integration moment versus vortex-integration moment such as that shown in fig- 
ure D5. The plot shown is for a wing swept 80° (8 = 10°) at an angle of attack of 30° (tan a/tan 6 = 3.27). 

Figure D5 shows that the wing-integration moments vary approximately linearly with the vortex- 
integration moments at the lower roll angles; this linear relationship breaks down only for absolute values 
of c|) greater than about 30°. The slope of the linear fit to the data between roll angles of -30° and +30° 
was therefore taken to be the value of 7 that the vortex-integration moments were multiplied by to get the 
hysteretic wing rolling moments. 
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(a) p = 5°. 


(b) p= 10°. 



a, deg 

(c) p = 15°. 


Figure D2. Variation of wing rolling moment coefficient with angle of attack for various wing sideslip angles. 







O A = 0.71 wing 



(a) Cj versus (3. 



(b) Q versus (j). 

Figure D3. Variation of wing rolling moment coefficient with sideslip and with roll angle for angles of attack below 
the stall. 
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Figure D4. Variation of rolling moment coefficient with wing roll angle obtained by integrating around wing and 
around primary vortices at a = 30°. 



^vortices 


Figure D5. Rolling moment coefficient obtained by integrating around wing versus that obtained by integrating 
around primary vortices at a = 30°. 
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Values of 7 were determined over a range of tan a/tan 5, as shown in figure D 6 . A cubic fit to these 
data (eq. (41)) was used to specify the reduction factor for a given wing at a given angle of attack in the 
wing rock calculations. 



0 1 2 3 4 5 

tan a/tan 6 


Figure D6. Variation of rolling moment reduction factor 7 with tan a/tan 5 . 
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Table 1. Computed Vortex Positions as Function of Wing Roll Angle 
for Different Angles of Attack 

(a) a = 10° 


<l>, deg 

XI 

y i 

x 2 

J2 

-52.5 

0.2804 

MBBm 

0.1073 

■BBS 

-45.0 

0.3585 


0.1274 


-37.5 

0.4289 

1.0125 

0.1538 

-0.7258 

-30.0 

0.4532 

0.9065 

0.1819 

-0.6394 

-22.5 

0.4262 

0.8063 

0.2026 

-0.6036 

-15.0 

0.3767 

0.7294 

0.2284 

-0.6027 

-7.5 

0.3272 

0.6764 

0.2537 

-0.6170 

0 

0.2975 

0.6843 

0.2975 

-0.6843 

7.5 

0.2537 

0.6170 

0.3272 

-0.6764 

15.0 

0.2284 

0.6027 

0.3767 

-0.7294 

22.5 

0.2026 

0.6036 

0.4262 

-0.8063 

30.0 

0.1819 

0.6394 

0.4532 

-0.9065 

37.5 

0.1538 

0.7258 

0.4289 

-1.0125 

45.0 

0.1274 

0.8331 

0.3585 

-1.1027 

52.5 

0.107.3 

0.9208 

0.2804 

-1.1781 


(b)a= 15° 


<t>. deg 

x \ 

Tl 

x 2 

T2 

-52.5 

0.3897 


0.1487 

mssM 

-45.0 

0.4981 


0.1800 


-37.5 

0.5832 


0.2220 


-30.0 

0.5844 


0.2571 

-0.5449 

-22.5 

0.5334 


0.2805 

-0.5462 

-15.0 

0.4682 


0.3008 

-0.5745 

-7.5 

0.4083 

0.6960 

0.3259 

-0.6093 

0 

0.3608 

0.6480 

0.3608 

-0.6480 

7.5 

0.3259 

0.6093 

0.4083 

-0.6960 

15.0 

0.3008 

0.5745 

0.4682 

-0.7622 

22.5 

0.2805 

0.5462 

0.5334 

-0.8559 

30.0 

0.2571 

0.5449 

0.5844 

-0.9804 

37.5 

0.2220 

0.6147 

0.5832 

-1.1191 

45.0 

0.1800 

0.7591 

0.4981 

-1.2310 

52.5 

0.1487 

0.8845 

0.3897 

-1.3201 
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Table 1. Continued 


(c) a = 20° 


<l>, deg 

*1 

y l 

x 2 

V2 

-45.0 

mmm 

-1.3655 

0.2358 

-0.5673 

-37.5 

wBm 

-1.1962 

0.2896 

-0.4154 

-30.0 


-1.0135 

0.3148 

-0.4200 

-22.5 

mUm* 

-0.8674 

0.3275 

-0.4730 

-15.0 

0.5215 

-0.7612 

0.3413 

-0.5301 

-7.5 

0.4525 

-0.6863 

0.3641 

-0.5812 

0 

0.4002 

0.6303 

0.4002 

-0.6303 

7.5 

0.3641 

0.5812 

0.4525 

-0.6863 

15.0 

0.3413 

0.5301 

0.5215 

-0.7612 

22.5 

0.3275 

0.4730 

0.6021 

-0.8674 

30.0 

0.3148 

0.4200 

0.6805 

-1.0135 

37.5 

0.2896 

0.4154 

0.7306 

-1.1962 

45.0 

0.2358 

0.5673 

0.6845 

-1.3655 


(d) a = 23° 


<t>. deg 

xi 

VI 

x 2 

T2 

-52.5 

0.6009 

1.2532 

0.2052 

IBB 

-45.0 

0.7913 

1.4581 

0.2848 


-37.5 

0.7683 

1.2370 

0.3320 

BBoi 

-30.0 

0.7193 

1.0364 

0.3480 

-0.3723 

-22.5 

0.6331 

0.8813 

0.3545 

-0.4483 

-15.0 

0.5482 

0.7692 

0.3651 

-0.5172 

-7.5 

0.4767 

0.6897 

0.3867 

-0.5753 

0 

0.4230 

0.6292 

0.4230 

-0.6292 

7.5 

0.3867 

0.5753 

0.4767 

-0.6897 

15.0 

0.3651 

0.5172 

0.5482 

-0.7692 

22.5 

0.3545 

0.4483 

0.6331 

-0.8813 

30.0 

0.3480 

0.3723 

0.7193 

-1.0364 

37.5 

0.3320 

0.3220 

0.7863 

-1.2370 

45.0 

0.2848 

0.4061 

0.7913 

-1.4581 

52.5 

0.2052 

0.7316 

0.6009 

-1.2532 
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Table 1. Continued 


(e) a = 25° 


<l>, deg 

*1 

y i 

x 2 

J2 

-52.5 

0.6699 

1.6241 

0.2260 

-0.6638 

-45.0 

0.8451 

1.5085 

0.3208 

-0.2951 

-37.5 

0.8150 

1.2601 

0.3595 

-0.2684 

-30.0 

0.7412 

1.0506 

0.3690 

-0.3451 

-22.5 

0.6512 

0.8903 

0.3714 

-0.4344 

-15.0 

0.5642 

0.7747 

0.3799 

-0.5101 

-7.5 

0.4914 

0.6923 

0.4006 

-0.5723 

0 

0.4370 

0.6292 

0.4370 

-0.6292 

7.5 

0.4006 

0.5723 

0.4914 

-0.6923 

15.0 

0.3799 

0.5101 

0.5642 

-0.7747 

22.5 

0.3714 

0.4344 

0.6512 

-0.8903 

30.0 

0.3690 

0.3451 

0.7412 

-1.0506 

37.5 

0.3595 

0.2684 

0.8150 

-1.2601 

45.0 

0.3208 

0.2951 

0.8451 

-1.5085 

52.5 

0.2260 

0.6638 

0.6699 

-1.6241 


(f) a = 27° 


<t>, deg 

Xl 

V 1 

x 2 

J2 

-52.5 

0.7651 

1.7346 

0.2571 

-0.5446 

-45.0 

0.8815 

1.5460 

0.3552 

-0.2014 

-37.5 

0.8389 

1.2808 

0.3862 

-0.2215 

-30.0 

0.7608 

1.0641 

0.3892 

-0.3211 

-22.5 

0.6679 

0.8992 

0.3876 

-0.4222 

-15.0 

0.5791 

0.7803 

0.3941 

-0.5041 

-7.5 

0.5053 

0.6953 

0.4140 

-0.5700 

0 

0.4504 

0.6298 

0.4504 

-0.6298 

7.5 

0.4140 

0.5700 

0.5053 

-0.6953 

15.0 

0.3941 

0.5041 

0.5791 

-0.7803 

22.5 

0.3876 

0.4222 

0.6679 

-0.8992 

30.0 

0.3892 

0.3211 

0.7608 

-1.0641 

37.5 

0.3862 

0.2215 

0.8389 

-1.2808 

45.0 

0.3552 

0.2014 

0.8815 

-1.5460 

52.5 

0.2571 

0.5446 

0.7651 

-1.7346 
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Table 1. Concluded 


(g) a = 30° 


<|>, deg 

x\ 

y i 

x 2 

n 

-52.5 

0.9161 

mssm 

0.3327 

-0.2283 

-45.0 

0.9164 

iSSl 

0.4021 

-0.0925 

-37.5 

0.8688 

1.3081 

0.4247 

-0.1606 

-30.0 

0.7873 

1.0830 

0.4185 

-0.2887 

-22.5 

0.6907 

0.9118 

0.4105 

-0.4060 

-15.0 

0.5996 

0.7884 

0.4140 

-0.4963 

-7.5 

0.5245 

0.6997 

0.4328 

-0.5672 

0 

0.4690 

0.6307 

0.4690 

-0.6307 

7.5 

0.4328 

0.5672 

0.5245 

-0.6997 

15.0 

0.4140 

0.4963 

0.5996 

-0.7884 

22.5 

0.4105 

0.4060 

0.6907 

-0.9118 

30.0 

0.4185 

0.2887 

0.7873 

-1.0830 

37.5 

0.4247 

0.1606 

0.8688 

-1.3081 

45.0 

0.4021 

0.0925 

0.9164 

-1.5864 

52.5 

0.3327 

0.2283 

0.9161 

-1.9155 


(h) a = 40° 


<t>, deg 

x \ 

vi 

x 2 

yi 

-52.5 

0.9932 

2.0498 

0.4799 

0.1583 

-45.0 

0.9831 

1.6812 

0.5384 

0.1248 

-37.5 

0.9468 

1.3884 

0.5469 

-0.0152 

-30.0 

0.8629 

1.1428 

0.5127 

-0.2057 

-22.5 

0.7574 

0.9536 

0.4831 

-0.3659 

-15.0 

0.6610 

0.8172 

0.4766 

-0.4787 

-7.5 

0.5834 

0.7175 

0.4915 

-0.5636 

0 

0.5269 

0.6381 

0.5269 

-0.6381 

7.5 

0.4915 

0.5636 

0.5834 

-0.7175 

15.0 

0.4766 

0.4787 

0.6610 

-0.8172 

22.5 

0.4831 

0.3659 

0.7574 

-0.9536 

30.0 

0.5127 

0.2057 

0.8629 

-1.1428 

37.5 

0.5469 

0.0152 

0.9468 

-1.3884 

45.0 

0.5384 

-0.1248 

0.9831 

-1.6812 

52.5 

0.4799 

-0.1583 

i 

0.9932 

-2.0498 
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Table 2. Nondimensional Vortex Strengths 
for Selected Values of tan a/tan 6 


tan a/tan 8 

Q 

0 

0 

0.25 

0.3366 

0.50 

0.4885 

0.75 

0.5687 

1.00 

0.6076 

1.50 

0.6558 

2.00 

0.6901 

3.00 

0.7415 

3.30 

0.7534 

4.00 

0.7818 

5.00 

0.8159 
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